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A model of quantum measurement is proposed, which aims to describe statistical mechanical 
aspects of this phenomenon, starting from a purely Hamiltonian formulation. The macroscopic 
■ measurement apparatus is modeled as an ideal Bose gas, the order parameter of which, that is, the 

amplitude of the condensate, is the pointer variable. It is shown that properties of irreversibility 
and ergodicity breaking, which are inherent in the model apparatus, ensure the appearance of 
definite results of the measurement, and provide a dynamical realization of wave-function reduction 
or collapse. The measurement process takes place in two steps: First, the reduction of the state 
. of the tested system occurs over a time of order h/{TN^^*), where T is the temperature of the 

O ' apparatus, and N is the number of its degrees of freedom. This decoherence process is governed 

by the apparatus-system interaction. During the second step classical correlations are established 
between the apparatus and the tested system over the much longer time-scale of equilibration of 
the apparatus. The influence of the parameters of the model on non-ideality of the measurement is 
, discussed. Schrodinger kittens, EPR setups and information transfer are analyzed. 
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1. INTRODUCTION 

C . 

Understanding the specific features of quantum measurements has been a long-standing question . Discussions 
on this subject and its interpretations started in the early days of the quantum theory Physical and philosophical 

reflections on the problem and consideration of its different aspects became a source for deep conclusions about the 
quantum world and its classical counterpart [^jjij. Also the current activity in this field |6-11| clearly displays its 



> 

QQ , complexity and multifariousness, and witnesses that its final settling is not yet close. 

' The main purpose of this paper is to exhibit the complete solution of a model for a measurement process, which 

, will show that two main paradigms of modern statistical physics, irreversibility and ergodicity breaking in phase 
transitions, are crucially relevant for the quantum measurement problem. Together with some other conditions they 

. provide all necessary ingredients for the realization of ideal quantum measurements. No additional postulates should 

■ be posed, since standard quantum statistical physics completely suffices for the self-consistent explanation of this 

"Y"; • phenomenon. 

' We will start with a discussion on the problem of quantum measurement and the main steps which were made in 

^ its interpretation and understanding. 
I 

o 
o 



1. General measurement 

Let us recall the general requirements that a measurement should satisfy. The quantity to be measured in the 
system S under study is represented by a hermitean operator x with eigenvalues Xk- The so-called pointer variable 
X in the measuring apparatus A may take values Xk in correspondence with Xk- The observation of Xk provides 
statistical information about the state of the system. Moreover, some measurements termed as ideal |^] can be used 
as filters to prepare the system in a state where x takes a well-defined value Xk ■ 

More precisely, we denote by p the density operator of the system under study, by R that of the apparatus, and 
by TZ the global density operator. To include the possibility of the most general descriptions of the system and the 
apparatus, we analyze the situation in terms of density operators rather than pure states. For simplicity we will 
operate in this section with discrete spectra. At the initial time t — just before measurement, the system is in an 
unknown state p(0), the apparatus is in a fixed state R{0) and these states are uncorrelated, so that the complete 
state is described by: 
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7e(0) = p(0)(8)i?(0). 



(1.1) 



The uncorrelated character of this state simply reflects the fact that for t < the tested system and the apparatus 
were not interacting. An evolution operator T transforms this initial density operator into the overall density operator 
TZ{6) at the final time 6 of the measurement, which should have the form 



r[p(o) ® 7^(o)] = 7^(0) ~ ^Pfc7^fc. 



(1.2) 



Here T does not depend on the initial state of the system, which is arbitrary and unknown, though it may depend 
on the initial state of the apparatus R{Q). By ~ we mean "as precise as desired provided that the parameters of the 
apparatus are tuned suitably" . 



In a precise and unbiased measurement the final possible states TZk entering Eq. (1.2) should constitute an orthogonal 



set, 



tr(7^fc7^^) ~ Sk. 



kl, 



(1.3) 



so as to ensure that they represent exclusive events to which ordinary (non-quantum) probabilities can be assigned. 
They may be distinguished from one another by observing the pointer variable X , which takes in each TZk the value 
Xk with a negligible statistical fluctuation, namely 



tr(7^fcX) = Xk 



tr(7^feX2) ~ Xi 



(1.4) 



Hereafter trA and trg indicate traces in the subspaces of the apparatus and the system respectively, while tr is reserved 
for the complete trace. 



Each Xk occurs in Eq. (1.2) with a probability pk, which is determined by the initial state of the system in the 
form 



Pfc -trs(p(0)nfe), 



(1.5) 



where H^ = |a;fc)(a;fc| denotes the projection operator onto the eigenvalue Xk of x in the Hilbert space of the system. 
Although the states TZk may depend on the initial state p{0) of the particle, the quantities Xk do not depend on p{0) 
for precise and unbiased measurements, but are determined by the structure of the apparatus only. On the other 
hand, the probabilities pk are determined by the initial state p{0) only. 

The fact that the process which leads from TZ{0) to TZ{9) describes a measurement is reflected in a specific feature. 
Due to the exclusiveness property (1.3), Eq. ( |l.2| ) represents the occurrence of a classical random quantity fc, with 
probability pk- It expresses that after the observation has taken place, the overall system is left in a state TZk with 
probability pk- Indeed, the density operator TZ{9) describes a statistical ensemble of measurements all performed 
under similar conditions rather than an individual experiment. Since the variable X has a definite value in each state 
TZk, one can count the frequencies of the different XkS and th us r ecover the unknown probability distribution pk- 
It is the specific type of correl atio n exhibited for each k in Eq. (1^), between the properties (L4) pertaining to the 



apparatus and the expression (1^) for the probabilities pk in terms of p{Q), which allows us to gain information about 
the tested system. 

The crucial problem of quantum measurement is therefore to explain how the evol ution process of the coupled 
state TZ{t) from TZ{Q) to TZ{9) can produce a transformation T which ensures Eq. (1^). The subsequent process of 
observation then merely amounts to the selection of a single term of Eq. ( |l.2| ) characterized by the value Xk of the 
pointer. This last step is by no means different from the analogous process in the classical probability theory [^^, as 
was stressed recently by van Kampen and one of us. It should be stressed additionally that observation and selection 
refer to the specific type of measurement process which has been previously performed. 



2. Ideal measurement 



The above measurement scheme is rather general, and in particular describes situations where the system itself 
does not have a definite state after measurement (e.g., a photon is destroyed when detected by a photo- multiplier). 

A theoretically and practically important class of measurements are the ideal ones, which leave the tested system 
as weakly perturbed as possible after selection of the result. In an ideal measurement, the possible final states TZk 
factorize as 
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(1.6) 



where pk is expressed in terms of the initial state p{Q) of the system as 



Pk = — nfep(o)nfe, 
Pk 



(1.7) 



and where the possible final states of the apparatus, characterized by the value Xk of the pointer variable as in 



Eq. (1.4), are orthogonal as in Eq. (1.3), 



JkU 



(1.8) 



and are independent of the initial state 7?,(0). Notice that the condition (1.8) is stronger than that given by Eq. (|l.3|) 



Later we will connect it with robustness of the apparatus as an information-storing device. 

Thus, after observation of the apparatus, and sorting of a given outcome X^, the system is prepared in the state 
(1.7) by means of the so-called "reduction of the wave-packet" or "collapse of the wavefunction" . Apart from bringing 
the system into an eigenstate of the observable x associated with the eigenvalue xu, the projection (^T^) does not 
affect its other degrees of freedom. 

A complete measurement theory should provide conditions under which an apparatus interacting with the tested 
system brings it into one among the reduced states p^. 



3. The standard approach 

When they deal with quantum measurements, textbooks usually justify the above properties by relying on the 
conventional arguments initiated by Bohr [PJ^Jl^. Somewhat qualitative and incomplete, this general line of reasoning 
became known as the Copenhagen interpretation. Different, though closely connected versions were summarized by 
Rosenfeld fl^ and more recently by van Kampen |^ among others. The first precise discussion on the measurement 
problem was given by von Neumann who clarified the issues, but was led to consider the properties of quantum 



measurements, in particular the reduction (1.7) as a postulate^ which complements the standard principles of quantum 
mechanics. This additional postulate is, however, not needed as one can show, using consistency arguments, that the 
reduction of the wave packet can be derived from the standard principles and from natural properties attributed to 
measurements, in particular, repeatability 

Nevertheless, a complete understanding of a quantum measurement requires its analysis as a dynamical process. 



A crucial point to be explained is the nonexistence in the final state (1.2) of so-called Schrodinger cat terms with 
respect to the index k. For simplicity let us specialize on an ideal measurement. Starting from any initial state of 
the system, the final density operator TZ{6) should commute with both the measured observable x of the system and 
the pointer observable X of the apparatus. This property ensures that the result of the measurement process can be 
described in the language of classical probability theory. Otherwise, if H-iO) did include off-diagonal contributions in 
k, the standard interpretation of the measurement process would fail. 

The Copenhagen school of thought overcomes this difficulty by saying that any physically acceptable apparatus 
has to be a classical system ||l|,||. Thus it cannot exists in a state with superpositions. More refined Copenhagen- like 
approaches ||l^,|l^,H state that the apparatus is a macroscopic system, and therefore coherent superpositions may 
arise, but cannot be detected at least when measuring certain observables. This situation was illustrated by an 
exactly solvable model ||l^ (see also [|l^ in this context), where a class of observables was proposed, which are indeed 
non-sensitive to superpositions. Nevertheless, a small modification of this class allows to produce superpositions 

In the von Neumann- Wigner approach the states of the apparatus are pure, as i?(0) — \^){^\, and Rk — 
\^k){'^k\i and the evolution is characterized by the mapping 

|Vfc)l*>^I^OI*fe> (1-9) 

between initial and final states of the compound system. If this evolution is applied to a coherent initial state of the 
tested system such as 



(1.10) 



then one finally gets the following state 'R-{9) associated with p(0): 
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n{e)=Y,ak\^k)\'i>k) a^(V^fc,|(*fe'l, (1-11) 

k,k' 



whereas the desired state 'R-{9) in Eq. ( |1.2| ) involves only the diagonal elements of TZ(9): 

n{9) = J2 l"fcl'l^fc)|*fc)(V'fcl(*fc|. (1.12) 

k 

Partial traces of TZ and TZ over either the apparatus or the system subspace are equal, but we should explain why the 
off-diagonal terms of TZ are never seen. Indeed, one has for the partial density matrix of the particle 

p{9) = trA^(e) = tTAn{9) = J2 IV'fc)(V'fe|- (1-13) 

k 

However, the situation described by TZ{9) does not correspond to any measurement, since the apparatus and the tested 
system cannot be in definite states with definite probabilities. Actually, since with modern experimental equipments 
one is able to detect mesoscopic and even macroscopic superpositions (see e.g., jr7| , p^ where recent results are reported 
in the context of charge and flux macroscopic superpositions in Josephson junctions), the question remains open and 
can be formulated in the following form: What arc the concrete properties of a system, that make it usable as a 
measuring apparatus characterized by the collapsed state Eqs. (1.1-1.7) ? 



4. Irreversibility 



An important requirement for the realization of Eq. (1.1-1.7) is the existence of irreversibility. Whereas (1.11) 



is a pure state, we expect the pointer variable to take some well-defined value Xk with the classical probability pk , 
and this implies the final state to be the mixture ( 1.12| ). In other words, the von Neumann entropy Sv't^(JZ{9)) 



-tvTZ{9) hiTZ{9) of the state TZ{9) in Eq. (1.12) is positive, while Svj<s{TZ{9)) = 0. More generally the von Neumann 
entropy of (1.2) can be shown to be different from that of ( |l . 1| ) , due to the elimination of the coherent terms in the 
final state. This implies a loss of information about the coherence, a loss which is required to ensure the classical 
interpretation of the measurement and the reduction of the wave packet. 

A measurement process should therefore be analyzed in the same way as an irreversible process in quantum statistical 
mechanics, a second reason for using density operators. The size of the apparatus should be sufficiently large, so that 
irreversibility and relaxation emerge from the microscopic reversible evolution. Otherwise measurements could n ot be 



ideal. The fact that the condition ( |1.2| ) cannot be realized with a unitary transformation from any initial state (1.1) 
is nowadays well-established with different degrees of generalization |§,|l9|. To understand this in simple terms, let us 
write down the ideal measurement transformation for two different initial states ^'■^•'(O), s = 1, 2 of the system: 

r[p(^'(0) ® R{0)] = TZ^'\9) ~ ^pi'Vfc ^ Rk, (1.14) 

k 

and assume for simplicity that the spectrum of x is non-degenerate: 11^ = \xk){xk\- Then unitarity of T requires 
trs[p(i)(0)p(2)(0)] = ti[TZ^^'>{9)TZ^^'>{9)], and conditions (0, |l|) imply: 

trs p(i)(0)p(2)(0) = J2{xk\p^^\0)\xi){xi\p^^H0)\xk) ~ Y.'^Xk\p^'^mxk){xk\p^'\0)\xk) (1.15) 

k,l k 

and hence 

Y,{^k\p^'\0)\xi){xi\p^'\0)\xk)^0, (1.16) 

k^l 

which cannot be true for arbitrary /9^^^(0) and p^^-'(O). Thus, the unitarity of T has to be gotten rid of. Actually, 
the situation is the same as in any relaxation process: The overall system is isolated, and the evolution of TZ{t) is in 
principle governed by Hamiltonian dynamics, but on suitable time-scales irreversibility occurs owing to the presence 
of a large number of degrees of freedom, which act as an external bath. Statistical physics is needed to explain 
this behavior, in which microscopic reversible equations of motion result in macroscopic irreversible ones. Just the 
same approach was followed recently by two of us, to discover that the standard issue of Brownian motion leads to 
incompatibilities with thermodynamics in the regime of quantum entanglement [po| . 
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5. Ergodic and decoherence approaches to quantum measurements 



A first application of statistical physics in support of the Copenhagen interpretation was given by Daneri, Longier 
and Prosperi ||l3|,|l2| . Somewhat related (but not equivalent) approaches are reviewed in Ref . [0 . After considering the 
measuring apparatus and measured system together as an isolated system, one attributes the absence of macroscopic 
superpositions to inevitable statistical uncertainties, which are present in macroscopic bodies. Mathematically this is 
reflected in different kinds of ergodic assumptions, which are reasonably creditable for those systems. However, these 
approaches have several drawbacks, which, in particular, originate from the fact that they do not provide dynamical 
mechanisms for the realization of quantum measurements. Extensive criticism on them can be found in |^]. 

There is another, nowadays not less influential, school of thought which we shall follow. It attempts to handle the 
problem also involving certain arguments from statistical physics [ pl|j2^ ,|9|, |Tl| . In this decoherence approach the 
loss of coherence is viewed as a process established by an external environment, which is generally understood as a 
collection of uncontrollable and unobservable degrees of freedom. A decoherence process suppresses superpositions of 
some special states, which are determined by the interaction between the environment and the system. However, in 
an ideal quantum measurement, the coherence associated with both the observable x of the tested system, and the 
observable X of the apparatus should disappear independently of the concrete form of environment-system interaction, 
whereas the other coherences which exists in p(0) should remain present in the reduced states pk defined by Eq. ( |l.7| ). 
The type of decoherence occurring in quantum measurements is thus very special and we shall relate its features to 
the apparatus-system interaction rather than to an environment-system interaction. 



6. Requirements on quantum measurement models 

Keeping in mind both the successes and shortcomings of the various existing ideas, we tackle in this paper the 
quantum measurement problem by investigating a specific model. The coupled evolution of the system and apparatus 
is treated as a dynamical process of quantum statistical mechanics. By deriving an explicit solution we wish to show 
how the various features of an ideal measurement, expressed by Eqs. ( |l.l[ - |l.7| ), can emerge from the microscopic 
dynamics generated by the Hamiltonian of our model. 

When choosing the model we have been guided by various conditions that an apparatus should satisfy: 

1) It should have a degree of freedom X which may relax towards definite values Xk- 

2) It should be macroscopic so as to ensure an irreversible relaxation. 

3) This relaxation should be selectively triggered by the interaction of X with the variable x of the measured system. 

4) The various values Xk should a priori be equally probable, so as to avoid any bias produced by the apparatus. 
Thus, the various final states Rk, characterized by the value of X^, should have the same entropy. 

5) The apparatus should be a stable and robust information storing device, which implies that the states Rk are 
nearly in equilibrium and that after the measurement has been completed, X is a nearly conserved collective variable. 

These properties suggest to take for the apparatus a suitably chosen macroscopic system which is able to undergo a 
phase transition, with Xk as an order parameter. Indeed, a phase transition is a macroscopic process with robust and 
stable (or at least metastable) outcomes. Notice that the existence of an order parameter implies ergodicity breaking 
in contrast to the purely ergodic view at measurement jl3| . 

6) The measured quantity x should be coupled to the order parameter, and the apparatus should amplify this signal 
received during its interaction with the system. This is achieved by noting that the value of an order parameter can 
be controlled by an infinitesimally small source. The microscopic variable will play the role of such a source, which 
controls X but otherwise does not affect the apparatus. 

7) The relaxation of the order parameter is ensured by its coupling with other degrees of freedom of the apparatus, 
referred to as a thermal bath. It is this coupling which, together with the thermodynamical limit for the apparatus, 
will ensure the specific type of relaxation discussed above. 

We shall work out a model, as simple as possible, subject to the above requirements. The tested system is a 
one-dimensional particle, and the quantity to be measured is its position. The apparatus is a non-interacting Bose 
gas, which has an easily tractable phase transition. The variable X is the amplitude of the condensate. This situation 
will be shown to be generalizable for an arbitrary tested system and an arbitrary measured observable (Appendix A). 
Although specific and not realistic the model is thus suggestive for more general measurements. The possibility of 
tuning the parameters will help us to find the limit in which the measurement is ideal and to explore some imperfections 
of the measurement. 
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This paper is organized as follows. In section 2 we present the model and discuss its equations of motion. Limits 
which are especially relevant for the quantum measurement problem, as well as exact solutions of the equations of 
motion in the Schrodinger and Heisenberg pictures are considered in section 3. In section 4 we show that the present 
model realizes the conditions of ideal mcasiiremcnts discussed above. There wc also consider characteristic times of 
this realization and discuss imperfections which arise due to an incomplete thermodynamical limit for the apparatus. 
Our conclusions are presented in the last section. Several technical questions are considered in Appendices A, B, C 
and D. 



2. THE MODEL 



1. The apparatus and its bath 

As a model for our apparatus we will choose a system of N non-interacting bosons in a three-dimensional cubic 
box with volume V and periodic boundary conditions. Its Hamiltonian reads 

HA = Y,^i4(^i^ (2-1) 

i 

where a ■ , are the creation and annihilation operators of each single-boson state, and Ej is its energy, given in terms 
of its wave vector k and of the boson mass M by: 

= (2.2) 

Notice that Bq = 0. This apparatus is an open system, namely it interacts weakly with a large external environment. 
If there were no other interactions, the Bose gas would relax with time towards the Gibbs distribution: 

PA = I eM-PHA + PfjM), Z = tv eM-PHA + PfjM), (2.3) 

where T = is the temperature and fi is the chemical potential. Both of these quantities are imposed by the 
environment. Finally J\f is the number operator 



Af=J2aUi- (2-4) 



In a realistic apparatus what we call the "external" environment is actually constituted by a large number of degrees 
of freedom which are part of the apparatus itself. Here we treat it as a separate thermal bath, which can exchange 
energy and particles with the Bose gas. We can imagine this bath itself as a Bose gas which is much larger than the 
apparatus, so that its intensive variables remain fixed once and forever. The bath is characterized by its temperature 
T, its chemical potential /i, and by a quantum coupling which induces a relaxation time to the Bose gas. This 
situation corresponds to the grand canonical ensemble for the apparatus although the overall apparatus-bath system 
is isolated. Whereas the value of T is fixed by the overall energy, the value of /i is determined by the overall boson 
number. We shall focus below on a condensed gas, with |/x|/T small as 1/y/N. The tuning of /x is then achieved 
through a control of the overall particle number. 

As a system examined by means of the apparatus, we take a particle living in one- dimensional space, with mass m, 
an external potential V{x) and Hamiltonian 

Hs = ^ + Vix). (2.5) 
The measured quantity is the position x, which is coupled to the Bose gas through the interaction 

Hi = -gxX, X = (4 + ao), (2.6) 

where g is the coupling constant. The quantity X attached to the apparatus will be our pointer variable. It is 
macroscopically accessible through observation of the density of the Bose gas provided Bose condensation takes place. 
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which requires low temperature and sufficiently small In this case we expect the signal x to be amplified, because 
will be extensive. On the other hand, a tested system interacting with all levels of the apparatus would certainly 
be less interesting, since the threshold of the infiuence on the apparatus would be diminished. We have chosen the 
simple interaction ( |2.6| ) for theoretical purposes, although it is not realistic since it can change the number of bosons. 
Altogether the Hamiltonian of the apparatus and the system during the measurement reads: 

H = HA-nAf + Hs+ Hi, (2.7) 

to which we should add the interaction of the Bose gas with the bath, see Appendix D. Since the number of bosons 



may change through exchanges with the bath and the observed system, we have included for convenience in Eq. (2.7) 
the contribution from the chemical potential. Later on we will show how to generalize this situation to an arbitrary 
tested system and measured observable, keeping the same measurement apparatus. 

2. Bose condensation 

As well known the three dimensional ideal Bose gas in equilibrium undergoes a condensational phase transition at 
sufficiently large density or small chemical potential. Let us briefly recall this phenomenon, since in our setup this is 
a crucial property of the apparatus as an information-storing device. We consider here the Bose gas submitted to an 
external constant field source J, which later on will be identified with the term ^^h/2gx in Hi: 

HB='^e^ala^- J{ao+al). (2.8) 

i 

We work in the grand canonical ensemble. When there is no source term, it is known that despite different 
magnitudes of fluctuations in the condensed phase, the grand canonical and canonical ensembles are equivalent for 



the non- interacting Bose gas |2^. However, the source term in Eq. (2^) can have a macroscopic effect only in the 
grand canonical ensemble. It controls the density of the condensate, which can vary owing to possible exchange 
with the bath. In the canonical ensemble where the overall density is given, the density of the condensate would be 
practically insensitive to the presence of the source. In our present situation the density, and not only the phase of 
(ao), appears as an order parameter controlled by the field J. This property is specific for our model of non- interacting 
bosons, and would be invalid for a realistic Bose condensate with interaction. 

At the equilibrium state with temperature T and (negative) chemical potential fi one diagonalizes the gibbsian 
density matrix by shifting the lowest energy operators as 

ao = ao-—, aj = aj--- (2.9) 



This leads to 



(ao) - (2.10) 
M 

<«o«o) = ;^^ + ^- (2.11) 

The averages concerning excited states will not change, since the field acts only on the lowest mode. For the total 
density of particles one gets 

N 1 V-/ t N (2M)3/2 de^ 11 J2 

— = — > aJaj = / -J7 \ 3 \ ^, (2.12) 

V V ' ' An^h^ Jo ePi^-t')-! y c"^" - 1 Vfi^' ^ ' 

Here we went to the thermodynamical limit, which, if x and thus J are well defined, means to make the following 
change: 

, , fd^xd^k ^ , X. /■ d^fc / ^ ,.(2Af)3/2 ^ ^ 
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We also separated out the contribution coming from the lowest state e = 0. 

We wish to measure x, which, if x is well defined, means that we wish to find the value of J = ^fijlgx through a 
macroscopic observation of the Bose gas. This is feasible by observing the total density ( ^.12 ), provided that the last 



term in this equation is finite in the thermodynamical limit, which requires: — s- when N ^ oo. In this case the 
particle density 



splits into a non-condensed part, 

A^n _ (2Af)3/2 

and a condensed part 



^ ^ ^ + ^ (2.14) 
V V V ^ ' 



r dsV^ _ 0.165869 AP/^ 



N T P 

which is our order parameter. Since N^/V is a known function of temperature, we can deduce Nc/V from the total 
density, provided that this total density is significantly larger than the critical value ( 2.15| ): 



In order to use the Bose gas as a measurement device for J, we shall require Eq. ( 2.16 ) to be dominated by its last 
term, which means that 

(2.18) 



M-of^l. (2.19) 



and that 

Under such conditions, the scalar source term x, which will be replaced later on by the coordinate of the tested 
particle, can be deduced from the density of the Bose gas and the characteristics of its bath: 

^,^2^(N^r^Y ^^^^^^ 



The sign of x results through Eq. (2.10) from the phase of the condensate. Low temperatures will improve the 



efficiency of the measurement, since they provide small values for ( 2.15| ) and thus increase the ratio N/N^ - Th e 



chemical potential of the bath should be fixed at some small value when N ^ oo, so as to ensure Eq. ( 2.1§| , 2.19 ) 
and g should be such that as g^x^Ti/ is of order N . 

The excited states of the Bose gas will not play a direct role in our model of measurement, since they are not 
coupled to the tested particle. Nevertheless, they contribute, together with the bath, to exchange bosons with the 
condensate, the density of which can thus be controlled by the chemical potential as well as by the source J . 

Notice also that, although the various macroscopic states characterized by different values of the order parameter 
Nc/V can be distinguished, they appear on the same footing, because both their entropy and their energy are the same 
in the thermodynamical limit. Indeed the contribution of the condensate to the entropy is IniVc, its contribution to 
(Hb) is 2/xiVc where ^ T, so that both become negligible as oo. This was required to prevent the apparatus 
from having an intrinsic bias. 



3. Equations of motion 
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a. Dynamics of the apparatus in its bath 



Before we examine the equations of motion of the overall system including the tested particle, the apparatus and 
the bath, we will investigate in this subsection the situation without the tested particle. At some remote initial time 
t = to the apparatus was in an arbitrary non- equilibrium state. At that time it starts to interact with the bath, and 
the arising dynamics for the apparatus will be described by means of a weak-coupling quantum Langevin equation p^ ] 
(see Appendix D for the derivation of this equation from the Heisenberg equation associated with the apparatus-bath 
Hamiltonian) : 



-[i/A - fJ.Af,ai{t)] - "fai + bi{t) = -i(wi -I- a)aj - 70^ -I- ^/2J bi{t), 



(2.21) 



where 



?iUJi = Ei, 



ha 



-fi. 



(2.22) 



Notice that the chemical potential fj, is negative, whereas the parameter a is positive. The Langevin equations are 
written for the Heisenberg operators of the apparatus only, but the presence of the bath is reflected through a friction 
term —70^, and a random Gaussian force operator bi(t), which satisfies: 



Mt),bi{t')] = s.,sit - 1'), Mt),b,it')] = [bUt),bUt')] = 0, 

{btmit')) = d.kS{t - t')n^^, {h{t)bu{t')) = {b\{t)bl{t!)) = 0, 
1 



(2.23) 
(2.24) 

(2.25) 



The most important consequence of Eqs. (|] 2^, ^.24 12^) is that they ensure relaxation with the characteristic 
time 1/7 of the apparatus towards the Gibbs distribution (2_^) with the temperature T and chemical potential /i 
imposed by the bath. This can be seen from the following exact solution of Eq. ( ^.21 ): 



-y{t-to)-i{uJi+a){t-to) 



ai(io) 



ds e-'^«-'('^'+")'^ b,{t-s). 



(2.26) 



In particular, all possible moments (aT "(<)a™(t)) calculated with Eq. (2.26) for t — io ^ I/7 are identical to those 



obtained through the Gibbs distribution ( |2.3D . For example the average number of particles (aj(t)ai(t)) in the level i 
evolves in time according to 



{al{t)a,{t)) EE n,{t) = c-'^^'-'«^n.,{to) + (1 - e-2'^(*-*»))n°^ 



(2.27) 



which shows that ni{t) relaxes to its gibbsian stationary value at the characteristic time 1/(27). 

Since as shown by Eqs. (2.6, 2.7) the tested particle interacts only with the lowest level of the apparatus, the 
equations of motion for the excited levels with i > 1 will be always given by Eq. ( 2.26 ). Therefore, in the further 
discussion we will leave these excited levels aside. 



b. Equations of motion including the tested particle 

Let us now consider the situation with the tested system. The interaction between the apparatus and the tested 
system is switched on at the initial time t = 0. For t < the overall initial state factorizes as in Eq. ( 
p(0) is an arbitrary state of the tested particle and i?(0) is the Gibbs distribution of the apparatus given by Eq. (^.3|). 
Indeed due to the assumed condition t — to 3> I/7 the apparatus had enough time to relax starting from any initial 
state at i = to- From now on we shall drop the index in ao = a, bo = b. The Heisenberg-Langevin equation of 
motion for the lowest level of the apparatus reads for t > 0: 

a = i[i/,a(t)] - 70+^2^ b(t) = -iaa + — ^ g x{t) -70+ ^/2^ b{t). (2.28) 

V 2?i 

This equation is solved exactly as 
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i{t) = e-^*-'"*a(0) + — ^ / ds e"'^^-'"'' x(t - s) + j ds g-^^"'"^ b{t - s). (2.29) 



For ^ 1, and when x is constant, Eq. ( 2.29 ) expresses t hat t he Bose gas relaxes towards an equilibrium state, 
where the particle n umbe r in the condensate is given by Eq. ( |2.16| ) with J — y/h/2gx, as ensured by the second term 
in the r.h.s. of Eq. ( 2.29 ). The average number of particles n{t) = {a'^{t)a{t)) in the lowest state is evolving as 

2 2 

n{t) = e-27t„(o) + ^^^^^^^^^ (1 + e-2Tt _ 2e-T* cos at) + (1 - c-^-'')n'"i. (2.30) 

The second term in the r.h.s. of Eq. ( 2.30|) is the contribution supplied by the source, which shifts the condensate 
density. 

Although the evolution of the apparatus when there is no source leads to a well-defined equilibrium state where 
n(0) = is large but not extensive as F ^ oo, the small interaction with the tested particle is sufficient to change 

macroscopically n at times t ^ 1/(27) if l/^l = is sufficiently small. This means that the apparatus together with 
its bath constitute a system which is nearly non-ergodic when Bose condensation sets in. 

The Heisenberg dynamics of the particle reads: 

F = ^[Hs,Fit)] - ±L[j:it),Fmait)+aHt)) (2.31) 
for any operator F. We find in particular 

i = —p, (2.32) 

m 

p^-V'{x)+ f dsx{t~s)x{s)+r^{t), (2.33) 







ry(t)-r;o(t) + ?7i(i), (2.34) 



T^o{t) = y/n^ g / ds e-'''{b\t - s)e'"" + b{t - s)e~'""), (2.35) 



where 





mit) = 5e~^*(a^(0)e-' + a(0)e--*), (2.36) 

X{t) = g"^ e"^* sin at. (2.37) 

The interaction of the particle with the Bose gas produce a force, which has a random part. The randomness of 
the noise r]{t) arises from two independent reasons: the statistical (uncertain) character of the initial state of the 
apparatus, which gives the contribution r]i{t), and the random character of b, b^ , which occurs through 770 (0- Recall 
that at t = the apparatus was in equilibrium at temperature T and chemical potential /z, 

(a(O)at(O)) = 1 + -^3^' («^(0)«(0)> = —p^- (2-38) 
Since b{t), b^{t) are themselves gaussian, ri{t) will be gaussian as well, with the noise autocorrelation: 

K{t, t') - {v{ty,v{t')) = Ko{t, t') + Ki{t, t'), (2.39) 

Ko{t, t') - Mt);Vo{t')) - ^ cos[a(t - t')] coth^ (e"''!*-*'! - e-^(*+*')) , (2.40) 

K^{t, t') = {m{t);vi{t')) - ^ cos[a(t - t')] coth^ e-''(*+*'), (2.41) 
where we define for any operators A, B: 

{A;B) = ^{AB + BA). (2.42) 

Notice that K{t, t') is time-translation invariant, although its separate parts are not. It is seen that for t + 1' 3> I/7, 
which corresponds to the stationary apparatus, only KQ{t — t') persists. The most important effect of the bath on 
the dynamics of the tested particle is the appearance of a new characteristic correlation time I/7 in addition to the 
time-scale set up by the maximal frequency of the apparatus and to the universal quantum correlation time h/T. 
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c. Brownian motion of the tested particle 



In order to compare Eqs. (2.32 ~ 2.37) with the standard quantum Brownian motion approach, one integrates 
Eq. ( 2.33 ) by parts. This yields an equation, 



1 /■* 

p = -V'{x) / dsx{t 

m Jo 



which has the usual form in terms of a friction kernel 

_ J a cos at + 7 sin at 



m 



■r 



(2.43) 



(2.44) 



and of a noise ?y(t). Notice, however, that the friction kernel does not have a definite sign with t. The last term in 
Eq. ( ^.43 ) renormalizes the potential. As far as one is interested in the state of the tested particle itself, the full noise 
acting on it is r]{t). However, for the global state of the particle and the apparatus r]i(t) is a deterministic object, 
and only 770 remains as noise. 



d. Validity of weak-coupling quantum Langevin equations 



When substituting the above white-noise quantum Langevin equations for the actual interaction between the ap- 
paratus and the bath, a crucial fact was that their coupling is weak, so that the damping time I/7 is much larger 
than both the (maximal) dynamical characteristic time td of the apparatus, and the characteristic correlation time 
of the bath h/T. Under these conditions it is possible, as shown in Appendix D, to introduce an effective quantum 
noise b{t) with white spectrum satisfying Eqs. ( ^.23 -2.25), and get the Gibbs distribution as the result of relaxation. 
The time td is expressed as 



1 

a 



h 

R 



as seen from the free part of Eq. (2.28). Altogether the parameters of the bath should satisfy: 



a > 7 or T > 1^1 > ?i7. 



(2.45) 



(2.46) 



where we have taken into account the upper bound ( 2.18| ) on Notice that excited levels of the Bose gas has a 
lower dynamical time 



1 



(2.47) 



It is expected that in a non-ideal Bose gas these characteristic times will influence also the lowest mode, making its 
characteristic dynamical time lower. The weak-coupling Langevin equation is discussed with more details in Appendix 
D. 



3. DYNAMICS OF THE MEASUREMENT PROCESS 



1. Approximate conservation of the measured quantity. 



The above equations describe the joint evolution of the tested particle and the Bose gas in the bath which ensures 
relaxation. This evolution will describe a measurement of a;, if the Bose condensate registers the statistical distribution 
of X in the initial state /3(0) under the conditions specified in the introduction. In the present section we will take into 
account the following condition for ideality of the measurement. 

If x{t) did change with time during the process, the final state of the apparatus would be determined not only 
by the statistical distribution of x{Q) but rather by the whole history of its change. So to be sure that we indeed 
measure the quantity a;(0), we will require that the characteristic time for x{t) to change appreciably from its initial 
value is much larger than the relaxation time I/7 of the apparatus. This time is itself shorter than the duration 
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9 of the measurement, since for t = the apparatus should display definite results through well-defined stationary 
states. Notice the conservation of x{t) does not mean that the state of the particle is not changing. One only needs 
approximate conservation of the quantity x. For the considered model this is realized if over the time scale 9 one can 
neglect the contribution of the kinetic energy /{2m) to the change of the coordinate x. Obviously, this is realized 
for sufficiently large m. The more precise condition m ^ g^9'^ will be found in Appendix C, where we discuss the 
deviations from non-ideality which arise from the motion of the particle during the measurement. For other systems 
this conservation of the measured quantity may result from a commutation relation: [_ffs, x\ = 0, as stressed by Wigner 
(see also |^ in this context). However, the actual mechanism is basically irrelevant provided that the time-scale 
for the variation of x is much larger than the duration of the measurement. 

Notice that taking the limit m — > cx3 as we do in the bulk of this paper does not mean that we are dealing with 
a classical particle. In fact, the initial density matrix {x'\p{Q)\x") of the particle is arbitrary, which expresses that 
the quantity x is subject to a fully quantum probability distribution at th e initial time owing to the presence of 
off-diagonal elements. The large mass enters the equation of motion ( 2.32 ), and it only implies that the diagonal 
elements {x\p{t)\x) remain unchanged during the measurement. 



2. Wigner function and intermediate Wigner function 

a. Wigner function 

For solving the dynamics in the Schrodinger picture, we find it useful to employ the language of the Wigner function. 
The description through Wigner functions is, of course, completely equivalent to that through density matrices, but 
frequently, and in particular for quantum systems with linear dynamical evolution equations, this allows to obtain 
and understand results in more easy way. This is because Wigner functions allows to use the classical intuition as 
much as possible, which appears to be insightful for such systems. 

We shall focus on the degrees of freedom x, p of the tested particle and a, of the condensate, leaving aside the 
excited states of the Bose gas which relax independently. To introduce a Wigner function not only for the particle 
but also for the apparatus degrees of freedom, we represent the annihilation and creation operators in the canonical 
way: 

a=i(X + iP), a) = -^{X -\P), (3.1) 



2fi ^2fi 

X^y|(a + at), P = iy|(at-a), (3.2) 

[X,P]=in. (3.3) 

Recall that the connection between a density matrix {x'\p\x") and the corresponding Wigner function w{x,p) is given 
for each degree of freedom as 

+ f IpI^ -\)= j ^ e'^'^' ^(^'^')' (3-4) 

w{x, p) = 1 de e-'«P/'"' -f I - |) . (3.5) 
Notice that in the present paper the normalization of the Wigner function is chosen as 

wix,p) = 1, (3.6) 

since the integration with this weight corresponds to a trace. The Wigner function of the particle and the condensate 
together will be denoted by yV{X,P,x,p); those of the particle and condensate separately will be denoted as w{x,p) 
and W{X,P), respectively. Obviously, one has 

'^^ W{X,P,x,p) ^w{x,p), f ^^W{X,P,x,p) = W{X,P). (3.7) 



2nh ' ' ' ' J 2nh 

The Wigner function at time t in the Schrodinger picture can readily be represented in terms of the Heisenberg 
operators. For example, for the tested particle the corresponding formula reads 
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w{x,p;t) = (trp(O) w{x,p;t)), 

w(x,p;t] = f — exp[—iax — ibp + iax{t) + ibpit)], (3-8) 

J 47r^ 

where x(t), p{t) are the Heisenberg operators of the particle, p(0) is its initial state, and the average is taken with 
respect to the full noise ?? = r/o + ■ 



h. Intermediate Wigner function 

We shall also find it convenient to use yet another formulation, where the degrees of freedom of the tested particle 
are left in the matrix representation, and the Wigner transformation (3.5) is taken only for the apparatus degrees 
of freedom. We will call this object the intermediate Wigner function, and denote it as V(X, P, a;', x"), where x', x" 
denote the corresponding matrix elements in the a;-representation, 

V{X,P,x',x")^ j die-'^P/''{X + ^,x'\n\X-^,x")= J dp e'P("'-"")/^>V(X, P, ^ ,p). (3.9) 

Here again the excited states of the Bose gas are left aside. 



3. Exact solution of the equations of motion in the Heisenberg picture 



We consider a hea vy (m — > oo), free {V{x) — 0) particle. The equations of motion for the apparatus have already 
been solved as Eq. ( 2.2(: ) for the excited states and Eq. ( 2.2£ ) for the lowest state. The latter equation, in which 
x{t — s) can be replaced by a;(0), is written in terms of X, P and h — as 



X{t) = X{0) e"^* cosai + P(0) e"''* sin at 



xiO)g 



- [a — (7 sin at + a cos at) e 

7^ + 



+ / ds e--^^ [e-'°"b{t ~s) + e'^^'b^t - s)] 
Jo 



P{t) = P(0) e^'^' cos at - X{0) e"'^* sin at 



72 + q;2 



[7 + [a sin at — J cos at) e 



-7*1 



+iv^ / ds e"'^^ [e'^^b^t - s) - e-'^'^t - s)] 
Jo 



For the tested particle we can solve Eqs. ( [2.32 , 2.32) when the operator x(t) does not significantly change as 

x{t) = a;(0), 



Pit) ^ piO) + xiO)m + / du?7(M), 



where 



Jo Jo 



(72 + a'^)at + (72 - a2)e"''* sinai- 2a7(l - e~^* cosai) 



(72 + q;2)2 



and x(^)j vi^) s-re defined by Eqs. (2.34 - 2.37). Hereafter the following formulas will be used: 

— c^''* coswt] — 7e~'''* sincji 



, 7[1 — e cos wtl + we '^^smut 
dse ^ cosujs = — 
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7^ + LO^ 



dse 







7^ -f up 



(3.10) 
(3.11) 

s 

(3.12) 
(3.13) 

(3.14) 
(3.15) 
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4. Exact solution of the equations of motion in the Schrodinger picture 



Since the above dynamical equations are linear, there is a direct connection between the Heisenberg picture and the 
Schrodinger dynamics in terms of the overall Wigner function for the tested particle and the apparatus. The equation 
for the common Wigner function of the particle and the lowest mode has the form 

W{X,P,x,p;t) = I dXo dPo dpo ^{X,P,x,p;t \ X^, Po,x,pa-Q)W{Xo, Po\Q)w{x,p^-Q). 



(3.16) 

Here we denote by Xg, PqiPo the variables of the Wigner function at the initial time. The variable x of the Wigner 
func tion r emains unchanged from the initial to the final state, since the Heisenberg operator x{t) is conserved. In 
Eq. (|3T6| ) we used the fact that the initial Wigner function W is factorized into the partial Wigner function of the 
lowest level W{Q) and that of the tested particle w{0): 



W{X, P,x,p;0)^ W{X, P; Q)w{x, p; 0) = w(a;, p; 0) - exp 

A 



2A' 



-X' 



2A 



-P' 



(3.17) 



where we took into account the fact that at the initial time t = Q the apparatus already relaxed to the equilibrium 
Gibbs distribution under the influence of the bath. 

To find the transition kernel $, we notice that the Heisenberg operators X(t), P{t) and p{t) of the Bose gas and 
the tested particle given by Eqs. (3.10 - 3.13) are linear and involve only the gaussian noise b^h^ . We can then use 
the method of Eq. (3.8) for the whole system. Recall that the moments of the Wigner function coincide with the 
corresponding symmetrized operator moments. Hence, $ has the gaussian form 



\/det B 

<^[X,P,x,p;t I Xo,Po,x,po;0) ^ (^2n)^/2 '^'^P 
where {Li, L2, L^) is the following three-dimensional vector: 



1 ' \ 
"2 ^ BikLiLk j , 

i,k=l J 



(3.18) 



Li = X - (X) = X - Xa e"'^* cos at - Po e"''* sin at 

P- {P) = P- Po e"'^* cos at + Xq e"^* sin at - 

gXo 



xg 



7^ + a"^ 



[a — (7 sin at + a cos at) e 



-7*1 



xg 



7^ + 



L3 = p - (p) = p - Po - xC{t) 



7^ + 



[7 + (a sin ai — 7 cos at) e 



[j + {a sin at — 'f cos at) e '''*] 

gPo 



7^ + 



(3.19) 
(3.20) 

a — (7 sin at + a cos at) e~'^*] . 

(3.21) 



Here we have used Eqs. ( 3.10| , 3.11, 3.13), and the average (...) is taken by the noise operators b, b'^ directly and 
through rjQ. The quantity ^(t) was defined in Eq. (3.14). Notice that since we are interested here in the common state 
of the particle and the apparatus, the term connected with 771 which gives rise to the last two terms in Eq. (3.21) 
appears as deterministic. The gaussian quantum noise enters through B, which is a 3 x 3 symmetric matrix with the 
following elements: 



[S-i]n = {{X - (X))2> = [B-%2 = {{P (P))') = 2 [1 - ^~""] coth 
[B-%2 = {X-{X);P-{P))=0, 

[B-^],3^{X~{X);p-{p))^(x-(X); f ds fjo{s) 



ha 
2T' 



[B-%, = {P-{P);p-{p)) 



P-(P) 



gh / 7 — 27e cos at + ^t* 
Y 72 + ^2 

gh / a ~ 27e~'*'* sin at — Q:e^27t 
72 + a^ 



coth 
coth 



ha 
2T' 
ha 
2T' 



(3.22) 
(3.23) 
(3.24) 

(3.25) 



[B-%3^{ip-{p))^)= f f ds,dS2Koisi,S2) 

Jo Jo 



q^h , ha 



27i+ 1 



-2^t 



72 + a2 



47(7e ^* cos at — ae sin at — 7) 
(72 + ^2)2 



(3.26) 
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If we adopt the following notations 



m 



then we can write the matrix B ^ as 



and hence the matrix B will read 
A + C' 



B 



1 
AA 



-c 

— (7 
1 



C{t) = 



B-^ ^\ \ 




1 
a 



A 









Eq. ( 3.18 ) for <i> can altogether be written in the explicit form: 



23 



A + + ^2 





1 



-^=_exp^--(X,+L,, 



A(t) 



(3.27) 
(3.28) 

(3.29) 




( 1 



(3.30) 



2AA 



(La - C^i - <jL2)' 



(3.31) 



where A, Ci cr, A, Li, L2 and L3 are defined by Eqs. (3.19 - B.28) as functions of time 



5. Measurement of a variable with a discrete spectrum 

The results reported so far were obtained for the measurement of the coordinate x which has a continuous spectrum. 
It is of clear interest to indicate how the obtained results can be generalized for other specific situations, e.g., for 
measurement of spin. Here we provide a simple remark, which will set our results in a more general context. 

Let us assume that the interaction between the tested system and the apparatus is still given by Eq. ( |2.6| ) but 
now the tested system is completely arbitrary, and x in this equation refers to one of its observables. In particular, 
it can have a discrete spectrum. For simplicity we still neglect the self-Hamiltonian of the tested system. In this 
general case, the complete Wigncr function for the system and the lowest level is no longer defined. Nevertheless, the 
intermediate Wigner function of Eq. ( |3.9| ) is still perfectly defined. Recall that this function employs (X, P) variables 
for the lowest mode, but uses the matrix elements {x' ,x") in the eigenrepresentation of the measured quantity x. 
As we show in Appendix A the intermediate Wigner function corresponding to Eq. ( 3.161 ) adequately describes the 
general situation that we consider. Though the complete density matrix of the measured system and the lowest mode 
might also be used, the intermediate Wigner function V is a more convenient object to deal with. 

An illustrative example for the measurement of an observable with a discrete spectrum is the spin-boson Hamiltonian 
H: 

Hi = \ga,X, (3.32) 

where the measured observable ^cr^ is the z-component of spin for the tested system. 

The situation of a discrete spectrum measurement will also be encountered below for our original model in spite of 
the continuity of the coordinate x. Actually, we shall consider an initial density operator p(0) involving two distinct 
values of x only, a situation which does not differ much from a genuine discrete spectrum. 



4. IDEAL MEASUREMENT: POSTMEASUREMENT STATES 



In this section we will be interested by the postmeasurement situation. Let us first resume the conditions of ideality 
of the measurement that we have encountered above. The relaxation time of the apparatus should be small compared 
to the duration of the measurement: 
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iO-> 1. 

The coupling with the bath should however be small on the scale of the dynamical time {2 At): 



7 ^ a 



h 



(4.1) 
(4.2) 



On the other hand, if we denote by i a typical value of the coordinate to be measured, 

= tr(p(0)x2), (4.3) 

its coupling with the apparatus should produce a finite condensate density, which according to Eq. (2.20) is expressed 
as 



0{N). 



(4.4) 



The fact that the bath ensures Bose condensation, but that the condensate density remains dominated by the coupling 
with the particle, imposes the condition ( 2.18 ), that is 

1 



1»^» . 
T N 



To fix ideas we shall assume in the following that 

T 



= O 



(4.5) 



(4.6) 



with T finite in the thermodynamical limit. This will imply, from Eq. (4^), that the coupling constant g is finite. 
Finally, if p is the characteristic value of the particle momentum, the approximate conservation of x during the 
measurement means that 



P 



< X. 



(4.7) 



or m ^ as shown in Appendix C. We expect that under these conditions both the tested system and the Bose 
gas have reached at time a quasistationary state that we wish to study. 

By using the condition (4.1) we first recall that the excited states become decoupled and thermalized with the bath, 
since their Wigner function at time 6 deduced from Eq. (2.26) reads 



coth 



2T 



exp 



1 



h coth 



2T 



(4.8) 



The transition kernel $ between times and 9 can be obtained from Eqs. ( |3.19 - 3.28, 3.31). Under conditions (4.1, 
4.2, 4.5, 4.7) its various ingredients reduce to 



m = t; coth— ~ -, 
2 Zl a 



m = 0, a{e) = 



2gh0 



L,i0)=X 



gx 



xg 



-Po, 



L3 - CLi -aL2=p-po + -{P - Po) 

a 



xg 



a 



where we have used £,{9) — g^9/a as follows from Eq. (3.14). We note that Eqs. (4.5, |4.6D imply 



repectively, and that Eqs. ([4.l| , [4.4|) imply 

m » -T, 



(4.9) 
(4.10) 
(4.11) 
(4.12) 

(4.13) 

(4.14) 



while gx/a is of order V Nh. 

The postmeasurement state of the tested particle and the apparatus will be investigated below in three steps. First 
we will discuss the partial state of the apparatus and that of the particle. Later we turn to the global state. 
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1. Apparatus 



As seen, the apparatus itself behaves as a Bose gas, subject to a source field proportional to gx{Q). Owing to 
Eqs. (4.12, 4.14) the variance AA of in Eq. ( ^.31 ) is much larger than (Nh/x)'^. We can therefore readily trace 
out the tested particle from Eq. (3.16, 3.31) by integrating over both p and pq. The factor w{x,pq; 0) of Eq. dsl^ ) 
thus generates the probability density: 



dpo 
2'Kh 



w{x,pq;0) = (x\p{0)\x) 



(4.15) 



for the coordinate x in the initial state of the tested particle. The resulting expression of the Wigner function of the 
apparatus at time 6 has the expected form 



WiX,P;0) = / dx {x\piO)\x) W^{X,P), 



(4.16) 



1 (X -9^^ ^ p2 



(4.17) 



The Wigner function Wx{X, P), where A = T/a, describes the quantum Gibbs distribution of the apparatus at 
temperature T and chemical potential fi = —ha, with a classical source J = y/hj2gx. 

For each possible value of the coordinate x of the particle the apparatus is thus in an equilibrium state, with an 
order parameter proportional to x. In spite of the quantum nature of the variable x(0) which is governed by the 
initial density operator p(0), it acts on the apparatus as a classical random object. We can understand this classical 
feature by noting that because a;(0) = x{t) during the process, a{t) and a^t) commu te w ith it: [a{t),x{0)] = 0, 
[a^t),x{0)] = 0. Therefore, the situation is very similar to that considered in section 2.2, with J = gx{0)y/h/2. 
However, there is a subtle point, since the field J is now random. Its quantum randomness arises from the initial state 
of the tested particle, which in general is not an eigenstate of the coordinate operator. Since the off-diagonal part 
of p{0) disappears owing to the large size of X{9)A{9), Eq. (4.16) shows that the quantum nature of the randomness 
is suppressed. It is seen as well that the resulting classical randomness is quenched, which means that all extensive 
quantities have to be calculated for a fixed field and then averaged at the last step. 

Notice that for the continuous spectrum a certain difficulty may arise if (for example) the initial state of the particle 
is an eigenstate of the coordinate: {x'\p{0)\x") = 5{x' — xo)5{x" — xq), in which case the quantity (4.15) diverges. 
There are several standard ways to overcome this difficulty Q. The simplest one is to consider a gaussian packet 
{x\'ip) = (27re)~^/'*exp[— (x — xo)^/(4e)] instead of a precise eigenstate of the coordinate. This state has a normalizable 
Wigner function. 



w{x,p] 0) = 2exp 



{x - xoY 2e p2 



2e 



(4.18) 



and we can let e — > in the last stage of calculations. Further on we will always assume that this procedure is implied 
when necessary. This justifies the replacement of (4.15) by 5{x 



Since the apparatus is a non-ergodic system, Eq. (4.16) means that it will occupy with a priori probabilities 
{x\p{0)\x) a state which is determined by the initial value of the coordinate. If the initial state of the particle is an 
eigenstate of the coordinate operator the field is not random, and there is only one state to which the apparatus can 
relax. So in this case we have a definite prediction, as it should be. 

In each state ( 1.17| ) of the apparatus, the momentum P associated with the lowest mode fluctuates around the 
value exactly as in the gibbsia n sta te imposed by the bath, but its coordinate X is shifted by xgh/\p\. This shift is 
of order V Nfi according to Eq. (4^), and it produces a finite shift in the density 



V 



2hV 



{x^ + p^-h) 



hg^x^ 
2p?V 



(4.19) 



of the condensate. The statistical fluctuation of X, equal to VA, is small compared to {X) oc %/ NU owing to the ffi'st 
condition in Eq. (4.13), namely A <C Nfi. 

The second condition in Eq. (4.13), namely A 3> ?i, entails that the variables X and P behaves as classical random 
variables, which is a natural requirement for the pointer variable of an apparatus. Compared to the shift {X) their 
fluctuations are of relative order N~^/'^ if we choose A = hO{^/N). 
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a. Amplification and registration 



We have just seen that the interaction of the tested particle and the apparatus results in a macroscopic change in 
the condensate density: It fixes the expectation value {X) to xgh/\^\ within fluctuations which are small in relative 
value. This large effect is a consequence of the condition ( [4.4[ ). The coupling g is sufficiently large to produce a shift 
in the number N^^ which is of the same order as the total particle number. However, it is sufficiently small so that 
the contribution of the interaction Hamiltonian Hi to the energy of the apparatus is negligible. 

The amplification of the influence of the tested system on the apparatus is due here to the smallness of The 
bath, which imposes on the apparatus the condition (4.5), prepares it before the measurement in a state where the 
condensate density is not yet finite, but where the smallness of makes the apparatus very sensitive to a source 
coupled to X. By making successive macroscopic observations of the value of X, one can then find the statistics of 
a;(0) through (a:|p(0)|a;). There is one-to-one correspondence between {X) and a;(0), without bias because the various 
values of the order parameter {X) yield identical values for the energy as well as for the entropy. 

A peculiarity of the model comes from the fact that the amplification factor depends on the coupling constant 
g and on the chemical potential of the bath. These quantities, which are both small, need to be known to let us 
determine a;(0) through {X). 

Moreover, if we wish to register the result of a measurement, which is the value reached by X at the time 9 when 
the interaction g is switched off, we need to imagine that the exchange of bosons between the bath and the apparatus 
is also switched off at the same time 9. The overall density of bosons as well the condensate density thereafter remain 
fixed in the apparatus, which is in ca nonical equilibrium after the time 9. Another theoretical procedure to freeze the 
condensate density at the value ( 4.19 ) would consist in switching off the coupling 7 before the end of the measurement. 
The fate of the tested particle after the time 9 is considered in Appendix B. 



h. Robustness 



We have explained how our apparatus realizes amplification of weak signals. This is only half of the way towards a 
good information storing device, because we yet should see whether another important property which is robustness 
is satisfied. In other words, if under influence of a weak fleld the apparatus has relaxed to a definite state, then what 
is the probability that it will leave this state spontaneously? If this transition probability is small, and can be made 
as small as it is desired, then the property of robustness is present. 

Let us assume that the apparatus has been brough into a state with 



a 



(4.20) 



In this state the apparatus has Wigner function Wx and density matrix R. We wish to calculate the transition 
probability to another state R' associated with x' under the effect of some perturbation. If these states were pure, 
the transition probability would read as usual: 



Pr(x x') = tv{RR'). 

For mixed states we use the same formula in terms of the overlap: 

dX AP 



Pr(x ^ x') oc iv{RR') = J ^^^^ Wx{X,P)Wx'{X,P). 



To be normalized this expression should be divided by Pr(a; ^ : 

g'^ix 



Pt(x 



exp 



Using Eq. (4.17) one gets 



4Aa2 



(4.21) 



(4.22) 



(4.23) 



It is clear that above the phase transition point, when is finite, this transition probability is of order one, so that 
no robustness is present as was to be expected. 

Let us consider the situation below phase transition where A ^ T/a. We then have 



Pr(2; 



exp 



4T| 



(4.24) 
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According to Eq. (4.5) the exponent is of order —N\^\/T. With Eq. (4.6) this exponent behaves as — a/ZV in the 
thermodynamical Umit, provided x and x' differ by a quantity which remains finite as TV — > oo. The probabihty 
therefore vanishes, as it should. 

The fact that the overlap between states of the apparatus associated with different values of x is negligible also 
expresses that different positions of the pointer variable constitute exclusive events. 



c. Accuracy of measurement 

The robustness reflects stability of the apparatus with respect to external perturbations. Another quantity, the 
accuracy, characterizes the strength of the noise due to the initial uncertainty of the tested particle, and due to 
spontaneous thermal fluctuations induced by the bath. One can estimate the accuracy of the measurement by 
evaluating the following quantity 



{(X^)}av-{(X)na 



(4.25) 



where the average (...) is taken with respect to the state of the apparatus, and where {...}av denotes the average over 
the initial distribution of the particle. This is the signal-to-noise ratio, and E <C 1 corresponds to a good measurement 
of X. Having used Eqs. ( 2.10| , 2.11 , 4.5), one finally gets 



hg'^{x'^}a 



(4.26) 



In the region where the condensational phase transition exists this quantity is small as T/ {\fJ-\N) — 0{N^^/^), provided 
that the thermodynamical limit is taken and that {a;^}av is finite. The accuracy is thus governed by thermal noise. 

Note also that apart from the above uncertainty the derivation of the measured quantity x from the pointer variable 
X by means of Eq. (4.19) involves the ratio g/^. The accuracy of the measurement is, of course, spoiled if the coupling 
constant g and the chemical potential of the bath are not controlled with precision. 



2. Tested particle 



Let us now consider the partial state of the tested particle. At time we can find this state by tracing out the 
apparatus from Eq. ( 3.16| ) using the approximations (4.E - 4.12] ). In this calculation we first note that the memory 
about the initial value Xq is lost. The variable Pq enters through the last term of the exponent of <i>, which at the 
time 9 reads: 



2AA 



(i3 - CLi - <yL2f = - 



xg 



(4.27) 



Since the apparatus is nearly in equilibrium at both times and 6, P and Pq are of order = ^/TJ^. We can thus 
neglect the term depending on the apparatus in the bracket of Eq. (4.27), because 



2AA, 



;(P-P^f = 



Ag'^T^O o? a 



1 

4^ 



(4.28) 



is small. This means that the overall system forgets about the initial state of the apparatus for 7^ 1. We can 
therefore readily integrate over the initial state W{Xo, Pq; 0) of the apparatus, then over the variables X and P, which 
yields 



w{x,p;9)^ / dpo w{x,po;0) 



1 



exp 



1 / g^e 



Due to the large value of AA the exponential factor in Eq. (4.29) is nearly constant. Indeed, we have 

AAi2 2hg^S:^ T 



76* = 76I 0{N^'^), 



(4.29) 



(4.30) 
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where x is defined by Eq. (4.3), so that the only effect of this exponential is to produce a cut-off which ensures the 
normalization of w{x,p] 9). Otherwise, w{x^p] 9) is practically independent of p, and its dependence on x is the same 

as that of the probability density {x\p{Q)\x) as expected. 

The density matrix in the x-basis associated with the Wigner function (4.29) is given by 

{x'\p{9)\x") = {x'\p{Q)W')^^V 
The large value of AA ensures that it practically reduces to the diagonal part of p(0) in the basis x. 



AA , 

7;{X 



■.Jin 

2|/iP ' 



(4.31) 



a. Decoherence time 



The above expressions for 'w{9) or p{9) show that the decoherence of the state of the tested particle with respect 
to X has been achieved at the time 9. 

In order to understand how and when this decoherence takes place during the interaction process b e twee n times 
and 9, we return to t he eq uations of motion of the tested particle which are given by Eqs. (3.S, 3.12 , 3.13 ). Both 
terms rjo and ryi in Eq. ( 3.13|) should h ere be treated as noise, since we eliminate the apparatus. The fluctuations of 
this noise are given by Eqs. (2.40, 2.41). Altogether we get for the Wigner function w{x,p; t) of the tested particle: 



w{x,p\t) 
where we define d{t) by 
d{t) 



dpo w(a;,po;0) 



y/2TTd(t) 



exp 



2d{t) 



{p-po-x^{t)f 



f ( Ah At2 K{ti,t2) = I Ah At2 Koih,h) + ^coth 

Jo Jo JQ 



ha f 1 — 2e cos at 
2T 



-27* 



7^ -I- 



and /p /„* Ah At2 KQ{h,t2) is given by Eq. (|3.26|) . By using Eqs. (|3.24i [3^2^ ) we find 



d{t)^[B-%3 + -[B-'] 

1 



13- 



Having translated to density matrix with help of formulas (3^, 3^) one gets the evolution as 



{x'\p{t)\x") = {x'\pm^")^^V 



2n 



2{x'~x"f+'-^{x'^-x"^) 



(4.32) 



(4.33) 



(4.34) 



(4.35) 



It is seen that the diagonal elements of the density matrix (those for which x' — x") are not changed at all, whereas 
the off-diagonal ones are damped with rate d{t)/h^ . In other words, the density matrix of the tested particle tends 
to the mixture formed by eigenstates of the coordinate operator. Indeed, the decoherence factor d{t)^ which measures 
how the density matrix is squeezed in terms of x' — x" , increases from to oo as the time t goes from to oo. Its 
asymptotic forms at short and long times are derived from Eqs. ( 3.24 , 3.26 ) as 



d{t)-- 
d{9) 



- — r, jt<$:at<$: i, 

a 

76I, j9 > 1. 



For long times 7^ ^ 1 the exponent of Eq. ( 4.35 ) behaves as 



2h- 



-ix'-x"y 



p- 



IMI 



76I = 76I 0{N^^^), 



(4.36) 
(4.37) 

(4.38) 



and we recover as in Eq. (4.31) the strong damping of the off-diagonal terms, which allows us to obtain of the proper 
state of the tested particle after the measurement. 

An analogous investigation for short times 'yt ^ at ^ 1 can be carried out using Eq. ( p^ ): 



(4.39) 
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This quantity can be much larger than 1 even in the region of short times t for which at = \fJ.\t/h <C 1, provided 



N-'/' « ^ « ^ = 0(7Vi/2). 



(4.40) 



The off-diagonal terms of p(t) thus disappear at the very beginning of the interaction of the particle with the apparatus, 
after a delay of order N~^^^. The characteristic time 



(4.41) 



over which the reduction of the state of the tested particle takes place is thus much shorter than the duration 9 of the 
measurement, since 



1 h 

7 m 



(4.42) 



The first inequality in Eq. (4.42) ensures the relaxation of the apparatus. The second inequality ensures that the 
stationary state will be gibbsian. The third inequality indicates that the decoherence of the tested particle takes 



place on a time-scale much smaller that the dynamical time h/\^\ given by Eq. (2.45). Note also that h/T itself is 
an important characteristic time-scale in quantum statistical physics, which characterizes the relevance of quantum 



versus thermal effects. As shown in Eq. (4.41) the macroscopic size of the apparatus reduces these quantum effects 
by a factor N^^'^. 

The reduction time r, which in textbook discussions is either taken to be zero or identified with the measurement 
time 9 itself, is definitely different from the latter in our model. Another interesting aspect of this difference is that. 



at those short time-scales given by Eq. (4.41), the change in the apparatus is still neghgible. Indeed, the variables of 



the bath have not yet changed at those scales, since the characteristic time over which the equilibrium between the 
bath and the apparatus sets in is I/7. On the other hand, the energy associated with the particle has as well not 
changed, because the influence of the kinetic energy can be neglected, and the coordinate x{t) remains constant. In 
fact we shall study in the next subsection the change in momentum associated with the second, imaginary term in the 
exponent of Eq. ( 1.35 ). According to Eqs. ( 4.36 ) and ( 4.45| ) below, this term is negligible compared to the first one for 
at <C 1 sin ce th e ir ratio for such short times is of order h£^{t)/d{t) ~ at\p\/&T — at 0{N~^/'^). Altogether, this means 
that Eqs. ( 4.39 , 4.41 , 4.42 ) provide an example of a situation where the reduction of the state of the tested particle 
occurs long before achievement of the measurement process, and without any energy cost, although this reduction is 
a consequence of the interaction with the apparatus. 



h. Backreaction of the apparatus 



The post-measurement Wigner function ( 4.29|) of the tested particle involves a shift xg^9/a of the momentum, that 
we have not yet discussed. The density matrix (4.31) accordingly exhibits oscillations within the small region in x' — x" 



where {x'\p[t)\x") is significant. Let us evaluate the order of magnitude of the corresponding average momentum: 



hxg'^9 hg'^x^9 a xh 



- — 76I = - - 76* OiN) > - 0{N). 



7 x^ 



X 7 



(4.43) 



The fluctuations of p, of order vAA = ^Jdiff)^ are also large, because the final state is localized in the x-space. 



They are expressed by Eq. (4.30). Although large, the value of {p) would be ineffective if it were smaller than these 
fluctuations. However, the ratio of the shift to the fluctuations is given at the time 9 by 



^^ = ^^ (^] ,9» 



^ ) 70 0{N^/^), 



(4.44) 



a large number in the thermodynamical limit. 

The shift of (p) is t heref o re an important effect of the measurement. Let us look how this shift increases at short 
times. By using Eqs. ( 4.32 , 3.14 ) we find: 



xm 



g axt^ 



7t < at < 1. 



(4.45) 
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At the time r when decoherence is being achieved, the shift (4.45) is of order 



\T J X X 



(4.46) 



We see that the change in {p) begins to be significant long after decoherence has taken place. This is consistent with 
the fact anticipated at the end of section 4. 2 a, that very little is yet changed in the apparatus at the time r. 

The ratio (4.44) is of order 1 at a time ri such that [a;^(Ti)]^ = d{Ti), where the shift becomes comparable with the 
fluctuation. Using Eqs. (4.36, 4.45) one finds 



r < Ti < 



(4.47) 



Here as in r, the time-scale is given by fi/T ^ but now the thermodynamical limit produces an enhancement. 

The large shift of {p) can be attributed to the interaction process which takes place between the apparatus and 
the particle in the time interval ri, 9. During this period, the particle acts upon the lowest mode so as to drive it 
towards a state with a finite density of the condensate. In response this large effect produces a boost in the average 
momentum (p). The phenomenon can be traced back to the second term in the r.h.s. of Eq. (2.33) for p. The factor 



x{t) which enters this term describes a deterministic effect produced on the particle by the apparatus in contact with 
the bath. The increase of |(p)| is thus a cumulative effect of friction, which accompanies the rise of |(^)|. 

Altogether the interaction of the tested particle with the apparatus produces on this particle two effects. It first 
reduces the state, suppressing the off-diagonal terms in x during the time r. Later on, between the times ti and 6*, i t 
yields a large value to the average momentum without spreading the distribution in x, as seen in Eqs. (|4.29| , 4.32 ). 
This second effect is probably connected with specific features of our model, namely, the choice of the apparatus and 
of its order parameter, and the form of the interaction Hamiltonian between the apparatus and the tested particle. 
It is also related to the existence of the continuous spectrum for x, which allows the rapid oscillations exhibited by 
Eqs. (^3l| , p5|) . 

One may wonder whether the large value of {p) reached at the time 6 is compatible with our hypothesis that the 
Heisenberg operator x{t) has remained practically constant over the time interval 0, 6. We show in Appendix C that, 
contrary to p(t), the equation of motion for x{t) contains no systematic drift term arising from the coupling with 
the apparatus and hence with the bath; it involves only a noise term which does not affect much x{t). We were thus 
entitled to neglect the variations of x between the times and 0. 



c. Einstein-Podolsky-Rosen experiment and speed of quantum signals 



The above analysis allows us to discuss an experiment of the Einstein-Podolsky-Rosen (EPR) type. Let us suppose 
that the tested system consists of two particles denoted by 1 and 2. They do not interact for t > 0, but they did 
interact in the past, which is reflected in an entangled wave- function of the tested system at the initial time t = Q: 



p(0) = |V)(^|, ^Y.'^k\xu)\Vk) 



(4.48) 



where \xk) are eigenfunctions of the operator x for the first particle, and \yk) are arbitrary normalized, not necessarily 



orthogonal functions in the Hilbert space of the second particle. As indicated by Eq. (4.18) a small dispersion should 
be allowed for x so as to normalize \xk)- The measurement of the observable x is realized as above, namely the first 
particle couples through its operator x with the apparatus as expressed by Eq. ( p.(^ ). However, the second particle 
does not interact with the apparatus. Eqs. (4.35, 4.39) take place as above with the slight difference that {x'\p{t)\x") , 
{x'\p{Q)\x") are matrices in the Hilbert space of the second particle. In particular, the reduction (collapse) of the 



initial state occurs on the time-scale predicted by Eqs. (4.39, 4.41), and it now provides: 



P{r) 



k 



\ak\^ \xk)\yk) {yk\{xk\- 



(4.49) 



The remarkable feature of quantum mechanics is that although only one subsystem is involved in the measurement, 
the total state of the tested system is reduced. 

This general analysis can be illustrated by the standard example of two spins produced in the singlet state by the 
decay of an object with angular momentum 0. The initial density operator (4.48) is 
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{siS2\p{0)\s[s'2) = ^Ss,+s.,fi ((5,i,y/s2,4 - Ss,,s'Js2,s[)- (4.50) 



As above we perform measurement on the z-component of the first spin only. For large N Eq. (4.39) leads to reduction 
Si = s'l after a delay r. This automatically implies S2 = s'2- 

(siS2|p(t)|s'iS2> - ^ ^si+S2,0 Ss,,s[ S,^,s'^. (4.51) 

Let us give some quantitative estimate of the characteristic reduction time-scale r. For a temperature T ~ 1 K we 
obtain 

T = lO-^TV-i/S. (4.52) 

Estimating N ~ 10^^ for a macroscopic system, we get r ~ 10~^^s. For a distance of 1 m between spins this would 
lead to a speed of order lO^^m/s. Of course, this does not mean that there is an information transfer at this speed, 
but only a change in our knowledge through p of the quantum correlations of the two spins. Indeed, as we noticed at 
the end of section 4.2a the energy of the system is constant for times of order r. 



3. The common state of the tested particle and apparatus 



In the evaluation of the common Wi gncr function ( ^3.16 ) of the particle and the apparatus, we use for the transition 
kernel $ the approximate expressions ( |4.9| - 4.12 ) as above. We recall that in the limit 'fO ^ 1 the variance of p, 



d{0) = Xie)A{9) 



(4.53) 



is large as well as X{9) — T/a. We saw at the beginning of Section 4.2 that this implies a loss of memory about the 
initial state W{Xq, Po;0) of the apparatus. In Eq. (3.16) we thus integrate over the initial Wigner function (3.17) 
and find 



where 



W{X,P,x,p;9) = / dpo ■^{X,P,x,p;9\x,po)w{x,po;Q), 



(4.54) 



*(X,P,x,p;6l I x,po)= 
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2A V aJ 



2A 2d{9) 



p-po~m^ + -P 



a 



(4.55) 



For d{9) o o the expression ( 4.54 ) reduces to the product of the Wigner function Wx{X, P) of th e apparatus, defined 
by Eq. ( 4.17 ) for each possible value of x, and that w{x,p;9) of the particle given by Eq. ( 4.29 ). This factorization 
merely expresses both the reduction of the initial state w{x,p; 0) into w{x,p; 9) and the registration by the apparatus 
of the classical random variable x. 

The finite size of d{9) in E g. (|4 55| ) entails a non-ideality of the measurement that we now consider. We shall find 
it convenient to rewrite Eqs. ( 4.54| ), (4.55) in two alternative forms. First we may use as indicated in Section 3.2b the 
density matrix representation only for the tested particle, while keeping the Wigner representation for the apparatus 
(intermediate Wigner function): 



V(X,P,x',x";9) = {x'\p(9)\x") - cxp 

A 



-±(x-^^-^ 

2A V a 2 



""^'-i-p2-^i(.'-.")F 
2A ah 



(4.56) 



where the density matrix p(9) reduces to a nearly diagonal form in x as expressed by Eqs. (4.31, 4.35| ). And finally 
the same expression can be presented in the complete density matrix representation: 



{X',x'\n{9)\X",x")^{x'\p{9)\x") 



/27rA 



exp 



1 

2A 



X' + X'- 



g x' + x" 
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2h' 
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-ix'-x"))' 



(4.57) 
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It is seen that the last smah term g{x' — x")/a in the r.h.s, of Eqs. (4.56) or (4.57) quantifies the entanglement, that 
is, the degree of quantum correlations between the apparatus and the particle. Indeed, if in Eq. ( 4.57 ) we neglect this 
factor, we again find that the overall density matrix simply factorizes into the contributions studied above separately 
for the apparatus and the particle: 



{x',x'\nie)\x",x") ~ {x'\R,\x") {x'\pm^"), 

where we check, using A = — TiT/ fi and fi — — ha, that 



= -L= f dX' dX" \X'){X"\ exp 
V 27rA J 
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2A 



X' + X" 



gx 
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A (X' - X"f 



(4.58) 



(4.59) 



is the gibbsian density operator of the apparatus for the pointer variable x — [x' + x") /2. This factorization, 
together with the fact that \p{0)\ describes the reduced state of the tested particle, shows that the program set up 
in the introduction by Eqs. (1.1 - 1.8) is achieved, and that our model describes an ideal measurement provided the 
various conditions (4.£ - 4.12) are satisfied. 

Let us turn to a more detailed discussion of the off-diagonal terms in the overall density matrix (4.57). As we 
already discussed in section 3. 5, Eqs. (4.56, 4.57) are valid for the measurement of any hermitean operator S of the 
tested particle, in particular those with a discrete spectrum. When dealing with this case, should be directly 

substituted by the eigenbase of S. Then the analogue of Eq. (4.57) reads: 



7^(^?) =^K|s,)(s,|(^i?,(0)+^exp 



-{S, - Sfe) - — 2-(s, - Sk) 



{Si\p{^)\Sk)\St){sk\ ® Rik, 



2h ' ' 2h- 

where pi = {si\p{0)\si) is the initial distribution of the measured quantity S of the tested particle, and where 



(4.60) 



R^k = -iL= [ dX' dX" \X'){X"\ exp 
v27rA J 
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satisfies 



trAi?i, 



exp 



W 1 



2h^ 



i^AR-k 



-^(x'-X"-%.-Sk) 
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(4.61) 



(4.62) 



For i ^ k Rik has almost the same form as Rn = Ri, given for the continuous case by Eq. (4.59), but with slightly 
shifted off-diagonal matrix elements. Notice that this shift is due to entanglement between the apparatus and the 
particle. 



It is seen from Eq. (4.60) that the off-diagonal terms of TZ are strongly suppressed with the exponential factor 

did) 



exp 



2h- 



-{Si - SkY 



(4.63) 



We again find the expected features for an ideal measurement, without the difficulties of the continuous spectrum 
discussed in Eq. ( 4.18[ ) and section 4. 2b . 

In order to understand the above suppression in terms of observables, let us imagine that one calculates the average 
in the state (4.60) of some observable T with matrix elements {X'\!Fik\X") in the Hilbert space of the particle and 
the apparatus: 



{T) = trA(.F,,i?,) + ^ exp 
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2h 



i^l - si) 
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{Si\p{0)\sk) tTAiJ^kiRik)- 



(4.64) 



To ensure that the second sum in r.h.s. of Eq. (4.60) is non-negligible in the thermodynamic limit, one needs: 

■die) 
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(4.65) 



at least for one pair (i ^ k). Using 
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\tTAiJ^kiRik)\'^ < tVA^ik^ki trA 



and Eq. ( |4.62| ), one can write the condition ([4.65|) as 



2T 

iYA^^^k^kt > -r- r exp 



(4.66) 



(4.67) 



No general principle prohibits the existence of such an observable T which will satisfy Eq. (4.67). However, it 
is needless to mention that in the considered large N limit it would be quite pathological. So, under reasonable 
conditions one only has the diagonal term in (4.64). 



The same conclusions as we just drew from the large size of d{9) do hold in the continuous case given by Eq. (4.57), 
except for divergences associated with the continuous spectrum, and we can rewrite Eq. (4.57) in the form 

(4.68) 



TZ{9) = y p{x)\x){x\(E)R^{d), 

where p{x) = (a;|p(0)|a;), which exhibits the form required by the ideal measurement conditions. 
Let us finally discuss for illustration two examples. 



a. Transformation of an eigenstate 



If the initial state \xi) of the particle is an eigenstate of coordinate, one has for the initial density matrix and 
Wigner function, 



{x'\p{0)\x") = {x'\xi){xi\x"), w{x,p; 0) = S{x - x^). 



(4.69) 



More precise normalization according to Eq. ( 4.1q ) provides a small width e to the (5- function in w and multiplies it 

by 



'Sne exp 



The state of the apparatus and particle for ^ 1 will be 

W{X,P,x,p; e) = W^,{X,P) S{x-xi), 



(4.70) 



(4.71) 



where Wx-^{X, P), given by Eq. (4.17), is the final Wigner function of the apparatus. As expected the measurement 
does not change the state of the particle; it leaves the particle and the apparatus uncorrelated, apart from the 
registration of the value xi in the latter. 



b. Decay of initial superpositions: may Schrddinger kittens survive? 

Let the initial state of the particle be a superposition of two different eigenvectors \xi) and \x2) of the coordinate 
operator, which appear with amplitudes ipi, ip2- For simplicity we will take these amplitudes real. At the initial time 
one has: 



{x'\pmx") = ^rfk{x'\x,){xk\x" 



(4.72) 



i,k=l 



The corresponding Wigner function reads, within the regularization of Eq. (4.18) which a provides a normalization 
factor (|4J0|), 



w{x,p\ 0) = ^ vfSlx - Xi) + 2Lpiip25 ( X 



4=1 



Xi+ X2\ p{xi — X2) 



cos 



(4.73) 



The physical interpretation of this formula is obvious. The first two, incoherent terms refer to localized states at Xi 
and X2- The cross term, which describes coherence, is localized half-way between xi and X2] through its oscillations 
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it is associated with quantum interference, a fact which is clear when one notices that it makes the Wigner function 
alternatively positive and negative along the line y = {xi + X2)/'2. In other words, the initial state is highly non- 
classical. 



Using Eq. (4.55) one finds the common Wigner function of the apparatus and the particle as 



WiX,P,x,p■,e)^Y.^^ W-.(^'P) 6{x-x,) + WifiX,P,x,p;e), 



(4.74) 



where the contribution from the interference term of w{x,p] 0) after integration over po is 

die) 



Wii{X,P,x,p-e) = 2^1(^2 exp 



2h^ 



(X2 - XiY 



Wy{X,P) 5{x-x)cos 



{xi ~X2) ( g g'^9 

p H P - X 

n \ a a 



, (4.75) 



with X ~ (x\ + X2) /2. The Wigner function (4.74) is a sum of two contributions W^ni and Wif . The first one, Wm, is 
the expression ( |4.6^ ) describing an ideal measurement. It is a positive function and just consists as expected of the 
incoherent mixture of two measured values x\ and X2 with classical probabilities i~p\ and •^\. 

The interference term Wit is strongly suppressed due to the factor (^(6*), which according to Eq. (4.52) yields in 
Eq. (4.75) an exponent of order ^dN'^l'^ provided \x2 — x\ \ is not very small. The disappearance of the contribution 
of the interference term expresses that the initially existing Schrddinger cats are automatically suppressed, so that a 
classical interpretation can be given to the final result of the measurement. 

For the continuous spectrum there can be cases where \x2 — a^il ^ Since the initial superposition is then 



small, this situation can be called Schrddinger kitten. As seen from Eq. (4.75), the decay of such a state becomes less 
efficient as |a;2 — a;i| decreases. We may thus wonder whether Schrodinger kittens could partly survive in a non-ideal 
measurement process. We no te, h owever, that the values xi and X2 can be separated in a measurement only if the 
transition probability of Eq. ( 4.24 ) is negligible, which requires 



(X2 - Xi) > 1. 



4T 



Since the exponent in the damping factor of Eq. (4.75) which characterizes the decoherence 



(4.76) 



(4.77) 



is much larger by a factor of order jdN than the exponent ( 4.76| ) which characterizes the robustness of the measure- 
ment, even the weakest Schrodinger kittens disappear in any measurement process. Any kitten that can be detected 
by distinguishing from each other the two interfering values of x, has the same fate as a cat: it does not survive. 



5. SUMMARY AND CONCLUSIONS 



In the present paper we have studied a simple model in order to get better insight on the question of quantum 
measurement. Our purpose was to describe in full detail the dynamical process which takes place during the mea- 
surement, due to the coupling between the tested object and the apparatus. As usual in this problem, we rely on the 
fact that the apparatus is a macroscopic object, so as to ensure the rapid decoherence which is needed to explain the 
classical nature of the interpretation of a measurement. However, in our approach, this type of irreversible behavior 
is merged with the idea that the apparatus should be able to evolve indifferently towards different macroscopic states. 
The selection of the outcome should be controlled by a small interaction with the microscopic observed object. The 
evolution of the apparatus should therefore be non ergodic, which we realize by identifying the pointer variable with 
an order parameter in a phase transition. The interaction with the measured system behaves as a small source which 
drives the actual value of the order parameter. 

We wanted our model to fulfill the requirements on ideal measurements listed by the end of the introduction. We 
also wished to be able to produce, in a consistent framework and from the first principles, using standard methods 
of quantum statistical mechanics, a full solution for the equations of motion which describe the dynamical process 
of measurement. This led us to choose an extremely simple model. Our apparatus is a non-interacting Bose gas in 
contact with a particle and energy reservoir, which can undergo a Bose-Einstein condensation. The pointer variable is 
the condensate density, which is sensitive to a coupling of the lowest-energy level of the gas with the tested microscopic 
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system. To fix ideas we have chosen for this system a one- dimensional particle, the position of which is to be measured, 
but generalization to other systems is straightforward. When the interaction between the system and the apparatus 
is switched on, the condensate density relaxes to a value in one-to-one correspondence with the possible values of the 
position of the tested particle. We find that the randomness of this position is directly reflected by the statistics of 
the possible outcomes. The off-diagonal elements of the initial density matrix of the particle are suppressed by the 
process, and only classical probabilities enter the description of the correlations between the initial position of the 
particle and the pointer variable. 

The various parameters of the model can be tuned, so as to explore the validity of the approximations which ensure 
that the measurement is ideal. The following requirements, which were expressed mathematically at the beginning of 
Section 4, are needed for ideality. 

a) The apparatus should be macroscopic. This large size plays a double role. On the one hand it ensures though 
decoherence the appearance of definite results in the measurement process. This means that in its final state the 
overall system composed by the apparatus and the tested system may be found in different mutually exclusive states 
with probabilities given by the initial distribution of the measured quantity. On the other hand, the macroscopic 
number of condensed bosons ensures a robust and accurate registration of the measured observable. 

b) Before the interaction with the tested system starts to act, the apparatus should be prepared in a state which is 
extremely sensitive to this interaction with the tested microscopic system. When the coupling is switched on at some 
initial time, the initial state of the apparatus becomes thus unstable, and it relaxes to another state determined by 
the measured object. This was achieved by letting the initial number of condensed Bosons be already large (as Vn) 
but not yet macroscopic. 

c) The relaxation time of the apparatus should be larger than the dynamical time. 

d) The coupling constant g should be finite, so that the source term produces a macroscopic effect on the condensate 
although the interaction term in the Hamiltonian is not extensive. 

e) The duration of the measurement should be larger than the relaxation time of the apparatus. 

f) The statistical distribution of the measured quantity, here the position of the tested particle, should remain 
constant during the whole measurement process. 

The model explains the collapse of the state of the tested particle as an effect of its coupling to the lowest level of 
the Bose gas. The thermal bath acts only indirectly, through the apparatus. We can therefore understand why this 
decoherence process eliminates the off-diagonal elements of the density matrix in the x-basis and not in another basis. 
Indeed, the quantum noise due to the environment affects directly the pointer variable X of the apparatus, and the 
tested system feels it only because it is coupled to X through the term —g x X oi the Hamiltonian. The association 
of the decoherence with the measured observable is thus a natural outcome of the model. Thus in our model the 
decoherence is determined by the interaction between the tested system and the apparatus. This interaction is a 
tunable property and can be controlled. Let us notice that in the standard decoherence approach P,pT[| this process 
depends on the interaction between the apparatus and its environment which is a hardly controllable quantity by the 
very definition of (unobservable) environment. 

Remarkably, the collapse of the state of the measured system takes place at the very beginning of its interaction 
with the apparatus, over a time t — h/ (TN^^'^^ which has the usual features of a decoherence time, proportional to 
h/T and small in the thermodynamic limit of the apparatus. This time scale should be contrasted to the much larger 
time scale I/7 associated with the relaxation of the apparatus. Once the state of the tested particle is reduced, very 
little has yet been changed in the macroscopic apparatus. It still takes a long time, of order I/7, for the apparatus to 
reach its new equilibrium position determined by the system. 

The reduction takes place for both the states of the tested particle and the apparatus, which remain only coupled 
by classical correlations at the end of the measurement. We have estimated the order of magnitude of the corrections 
to this ideal situation, and seen that they become extremely small as N increases. In particular the broadening of 
the density matrix of the particle around its diagonal elements is of order less than N"^/"^. Related to this aspect 
is the suppression of Schrodinger cats (interference effects of states located at two different positions), and even of 
Schrodinger kittens (similar states at two nearby positions), which cannot survive a robust measurement. 

If the tested system involves degrees of freedom other than the one which is measured, our model shows that they 
remain unaffected by the process. As expected for an ideal measurement, the density matrix of the system is then 
simply projected according to Eq. (1.7) by the interaction with the apparatus. In particular, in Einstein-Podolsky- 
Rosen setups, our analysis confirms that the measurement of one particle implies the collapse of the full system, thus 
also of constituents that are spatially separated. This means that quantum information is transferred at a large speed, 
which depends on the size of the apparatus, and can thus be extremely large. 

Although oversimplified, our model has many generic features which are expected to occur in realistic measurement 
processes. However, taking advantage of the lack of interactions in the Bose gas, we have used an order parameter 
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which is the density of the condensate. This is a pecuHar property, since in a real Bose gas only the phase of the 
condensate, not its amplitude is an order parameter. A drawback of this situation is the fact that the back reaction of 
the apparatus on the system is stronger here than expected for a truly ideal measurement: Indeed we have seen that, 
apart from the reduction of the state to a nearly diagonal density matrix, a large amount of momentum is transferred 
by the apparatus to the particle. Another drawback of our model is that it simulates, rather than completely involves 
non-ergodicity: Although very sensitive to perturbations, our initial state is stable. It would be desirable to work 
out a more elaborate model where the initial state is metastable, and can be displaced towards several possible truly 
equilibrium states characterized by an order parameter, this displacement being controlled by the interaction with 
the system. 

Other difficulties have been encountered above due to the continuity of the spectrum of the measured quantity. For 
instance, the large back reaction on the particle momentum is related to this continuity. However, such difficulties 
are not an artifact of our model as continuous spectra are known to cause difficulties in many other circumstances. 

APPENDIX A: SOLUTION OF EQUATIONS OF MOTION FOR AN ARBITRARY MEASURED 

QUANTITY 



Here we discuss the solution of equations of motion ( 2.2g , 2.31 ) for an interaction Hamiltonian 



Hs - ~gXS, (A.l) 
where S is an arbitrary hermitean operator, which can in particular have discrete spectrum. In the present appendix 



we will distinguish operators by a hat. Let us rewrite Eq. (2.28) in terms of X and P: 



^^X = aP-^X + fx{t), fx{t) = ^(P{t) + h{t)), (A.2) 

^P = -aX^^P + gS + fp{t), fp{t)^x^(b\t)~b{t)). (A.3) 
at 

For a general operator of S no Wigner function exists in its standard sense, but for us it will be enough to operate 



with the intermediate Wigner function V(X, P, s',s") defined in section 3.2 b. Recall that this object is a Wigner 
function with respect to the lowest level of the apparatus (here the corresponding Wigner function can be defined) , 
but is a density matrix in the eigenrepresentation {|s)} of S'. 

To write the equations of motion in terms of V we will use a method described in |27j , which replaces the operator 
equations (A.2), (A.3) by the stochastic equation: 

dtV{X, P, s', s") = -dx{[aP - -iX]V) + dp{[aX + 7P]V) 
-dx{fxV)~dp{fpV) 

-^^^9^pV + \x{s'~s")V. (A.4) 

Here fx{t) and fp{t) behave as classical noises, which have exactly the same average and autocorrelation as the 
corresponding quantum quantities fx{t) and fp{t) after the symmetrization of Eq. (2.42). The true intermediate 
Wigner function V is obtained from V by averaging with respect to these noises fx [t) and fp {t) , 

V(X,P,s',,s") = (V(X,P,s',s")) (A.5) 

The first line in the r.h.s. of Eq. ( |A.4| ) is the standard drift contribution of the Liouville- Wigner equatio n. T he last 
line in this equation is as well explained rather simply: This is just the result of the Wigner transformation ( |3.4[ ) which 
was applied for the apparatus to the corresponding term {i/h){s'\\H^ p\ \s") in the density matrix representation. 



Now it is easy to see by direct substitution that the solution of Eq. (|A.4| ) reads: 

ft 



V{X,P,s',s") = exp 



|(/-.") / dt'X(t') 



S{X-X{t))5{P-P{t)), (A.6) 



where X(t) and P{t) are the solution of the following c- number equations 
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fx it), 
s' + 

-aX -jP + g— 
This means that one has the following solution for V: 



At ' 



V(X,P,s',s") 



d^od-Po ( exp 



15, 



{s'-s") / dt'X{t') 



fp{t). 



(A.7) 
(A.8) 



S{X - X{tmP - P{t))j W{Xo, Po){s'\pms"), (A.9) 



where W{Xq, Pq) is the initial Wigne r function of the lowest level, and Xq = -'^(O), Pq = -P(O) are initial values of 
X{t), P{t) which are inherent in Eqs. (A/7, A^). Since the exact solution of thes e equ ations is available in the form of 
Eqs. ( 3.10| , |3.11 ), a little patience is suffi cient t o ver ify that V{X, P, s' , s") in Eq. ( |A.9[) coincides with the intermediate 
Wigner function corresponding to Eqs. ( 3.16 , 3.18 ) provided that one makes the identification s' ^ x' , s" x" . 



APPENDIX B: THE STATE OF THE PARTICLE AFTER MEASUREMENT 

Here we discuss the post-measurement evolution of the state of the particle. We will assume that at time t ~ B the 
interaction between the particle and the apparatus has been instantaneously switched off. The switching is needed to 
ensure that the measurement will remain ideal, whereas its instantaneous character is taken for simplicity. Indeed, if 
the tested particle is interacting with the apparatus long enough, its coordinate will start to change due to its own 
Hamiltonian iJsi as well as due to interaction with the apparatus. Since the apparatus is itself interacting with the 
bath, sufficiently long interaction of the particle and apparatus will finally lead to relaxation of the particle towards 
certain steady state, which is independent of its initial state. This will violate the condition of ideality. 

Therefore, the time d was assumed to be much larger than I/7, so that the apparatus has enough time to relax to 
its stationary state and monitor the results of measurement. On the other hand, 9 was assumed to be small enough 
so that effects connected with change of x are not yet relevant. For t > 6 the tested particle thus follows its free 
evolution, which is described by free Heisenberg equations: 



x{t) ^ x{9) + 



-p{9) 



(B.l) 
(B.2) 



Due to the instantaneous character of the switching, p{9), x{9) arc those operator values which the particle reached 
during the interaction with the apparatus. 

The dynamics from i = to t = 6* is described by Eq. ( 4.32| ), and ioY t > 9 one has 



w{x,p;t) 



dxo dpo w{xo,Pq,9)S{p - po) S [ x - xq 



-Po 



(B.3) 



Let us now investigate w{x,p; t) at t > 9 for the initial state at t = given by Eq. (4.73): 
2 



){x,p\t) ^^Lpfsix-x. 



k=l 



t- 



-p 



, t~9 X1+X2 

+2tfiip20 X p 



exp 



d{9){xi^X2y 
2h^ 



2% 



(p-m^) 



(B.4) 



It is seen that the last contribution to this equation is due to incomplete reduction. If it is not suppressed totally in 
the course of the measurement (because N is not sufficiently large) , it will persist during further evolution, and might 
in principle be observed. 



APPENDIX C: MOTION OF THE PARTICLE DURING MEASUREMENT 



The purpose of the present Appendix is to discuss what happens if the change with time of the measured quantity 
X is not negligible, and if it cannot be treated as a constant of motion. In practice the conservation of the measured 
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quantity is ensured only over short times. However the duration of the interaction with the apparatus should be 
sufficient so as to ensure registration, and there will arise a source of non-ideality. Our purpose here is not to develop 
a full account of this non-ideality, but just to display on which time-scales its presence is not relevant. 

We will investigate Eqs. (2. 28, ^.32 , ^.32 ) on times where changes of the measured quantity x become noticeable. 
To keep the situation free of instabilities, we will make a natural assumption that the tested particle is subjected to 
a confining potential. For simplicity this potential will be taken to be harmonic: V(x) — mui'^x^ /2. 



Dynamics 



The general solution of Eqs. (2.32, 2.33) is obtained with help of Laplace transformation. Recall the following 
standard relations between functions A{t), B{t) and their Laplace-transforms A{s) — dt e^'**A(f) and B{s) denoted 
in this appendix with a hat: 



C 



1^ dt'A{t - t')B{t')^ = A{s)B{s), C |i} = ~A{0) + si(s), 



where A = f^ A. Thus, the solution of Eqs. (|2.32| , |2.33D reads 



= -fisMs) + im + sx{0))f{s), 



where 



is the Laplace transform of (sec Eq. (2.37)). Finally one has 



x(t) = x{Q)f{t) + -p{0)f{t) 



1 



dt'fit - t')7j{t'), 



p{t) = p{0)f{t) + mx{0)f{t) + / dt'fit - t'Ut'), 

Jo 



(C.l) 

(C.2) 
(C.3) 

(C.4) 

(C.5) 
(C.6) 



where f{t) is the inverse Laplace transform of f{s). 

Substituting Eq. ( p.5| ) into Eq. ( 2.29| ), we obtain the solution of the Heisenberg equation for the pointer variable: 



X{t) = <?x(0) / dt'e-'^'' smat'fit - t') + / dt'e"^'' sinat'f{t - t') + -C'^ {(—f + l] 

Jo m J„ 9 l\m J } 



(C.7) 



Its last term, where rj is given by Eqs. ( 2.34 -2.36) describe s the effect of the noise b, 6^ and the remanence of the initial 
conditions X{0),P{0) of the apparatus. We recover Eq.( 3.1C ) in the large m limit, for which /(t) = t if uiQt <C 1. 
However, instead of being controlled by x(0) only, the order parameter X(t) now depends on the initial state of the 
particle through both a;(0) and p{0). Since due to the uncertainty relation {p{0)'^) (x{0)'^) > Ti^/A these quantities 
cannot have definite values simultaneously, {X{t)) is always fluctuating with the initial state of the particle. This 
violates conditions ( |l.4 P^ ) according to which if the tested system starts its evolution from one o f the eigenstates 
of the measured quantity, then the result displayed by the apparatus can be definite. Moreover Eq. (C.7) shows that 
the statistics of X(t) is governed by the statistics of a;(0),p(0), or equivalently by the full density operator p(0) of the 
particle at the beginning of the measurement, including off-diagonal elements. It is the disappearance of the term in 
p(0) in Eq.(3.10) which allows X(0) at the end of the measurement to depend only on the diagonal elements {x\p{Q)\x). 
As we shall see below, this occurs for m ^ 9^0^ , that is, for a sufficiently short duration of the measurement. Once 
more this shows that the short-time limit is a necessary condition for (nearly) ideal measurement of the coordinate. 

Remember that the Heisenberg equation (2.28) for X{t) depended only on the position operator x{t)^ not on the 
momentum operator p{t). Its solution given by Eq. (2.29) expresses the action of the particle on the apparatus as 
a memory effect depending on x{t) at earlier times. Using the equation of motion (C.5) for a;(i), we have here re- 
expressed X{t) in terms of a::(0) and p(0). The momentum of the particle thus came out from the elimination of the 
history of x{t). 
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Short-time expansion 



To investigate the short-time behavior of the tested particle in more detail, we adopt the following large-mass 
approximation: 



fis) 
fit) 



1 



, 1 m 

+ lJq (s^ + o^q)^ to 



3 , di'[sin(tJot') - t' ujq cos ujot'] sm{a[t - t'])e"''(*"*'^ 



where 



(C.8) 
(C.9) 

(C.IO) 



is the characteristic frequency connected with the mass of the particle and the interaction with the apparatus. This 
frequency will be assumed to be the smallest characteristic frequency in our problem. Eq. ( C.9 ) is a short-time 
expansion, which is valid for 



(C.ll) 



Expansion ( |C.9D , when substituted into Eq. (C.5), produces the coordinate as a sum of two terms: The first one is 
due to the free motion in the potential V{x), and the second one represents a deterministic correction arising from 
interaction with the apparatus. The interaction with the apparatus should be switched off before this term becomes 
comparable with the first term in Eq. (|C.9D. In particular, if t is so small as t ^I/ivq we obtain: 



m = t+ 



Jo 



])e- 



-i(t-t') 



(C.12) 



When this equation is substituted into Eqs. {G.t, |C.6|), it is seen that in all terms besides the second term in the 
r.h.s. of Eq. (C.6) the correction to f{t) = t in E q. (|C.12 ) can be dropped under the condition ( p. 11 ). If one will 
take additionally: x{Q) > p(0)i/™, Eqs. ^ p.l2D produce Eq. ( p^ . 

It is important to notice that in this approximation there is a deterministic influence of the bath on the momentum. 
This is just friction, which may enhance or reduce (p). However, there is no such a systematic influence on the 
coordinate. It was therefore legitimate to assume in the bulk of this paper that in the Heisenberg representation we 
have xit) ~ a;(0) in the time interval 0, 9 in spite of the large value reached by {p) at the time due to friction (Section 
4.2b). 



Relaxation of the tested particle 



Although Eq. (C.9) has been sufficient for our purposes, we will mention here how the full f{t) behaves. This will 
allow us to understand the long-time behavior of the tested particle. One has to obtain the roots 71, 72, 73, 74 of the 
following equation 

+ ilk - lf][ll + - ^ (C.13) 

which provide 

4 

/(t) ^Y^AkC-^"', Ak = Lim,^_„,(s + 7fe)/(s). (C.14) 
fe=i 

Using the smallness of $7 the above roots can be obtained approximately as 

71 7 + m -I- (5i, 72=7i, (C.15) 
73 = iujQ + (52, 74 = 73' (C.16) 
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where, after using a ;:!> 7, a ^ luq, we find: 



^ 1 

2i (7 + ia)2 + ^2 

an^ 1 



2iwo (1(^0 — 7)^ + oi' 



2^' 



(C.17) 
(C.18) 



As seen from Eqs. ( C.15 , C.16 ) the particle has two relaxation times: I/7 (which is also the relaxation time of the 
apparatus) and a much longer one l/(Re ^2) = a^/fl^j. They are widely separated since a^/fi^ 3> 1. After the time 
/n^-f = c?g^ jm^ all information about the initial state will be forgotten by the particle if it still interacts with the 
apparatus, and it will relax to its equilibrium state imposed by the bath temperature. 



APPENDIX D: WEAK-COUPLING LANGEVIN EQUATION 

In the present appendix we will derive from a consistent Hamiltonian formulation the Langevin equation (2.21) 
for the apparatus-bath dynamics. More information on weakly coupled (weakly damped) dissipative systems can be 
found in [24]. 

For simplicity we deal only with the lowest energy level of the apparatus, the Hamiltonian of which reduces to the 
chemical potential term —[ia)a = ?iaa^a; as in Eq. (2.28) we drop the index 0. The dynamics of the excited state 
operators would be obtained by substituting a + 0;^ to a. The considered mode is coupled to the thermal bath, 
which is also a non-interacting Bose gas having single-particle levels to with creation operators i^a ^^'^ energies h£lrn 
(including the chemical potential of the bath) . This bath is equivalent to the dense set of harmonic oscillators usually 
considered in quantum Brownian motion [20, 24, 25], but here the coupling should account for transfers of bosons 
between the bath and the apparatus proper, with conservation of the total particle number. We can thus describe 
the apparatus with its bath by the total Hamiltonian 

The resulting equations of motion in the Heisenberg picture are: 

a = -iaa - i^ CmCm, (D.2) 



irn = -if^mCm " icJ^^O. (D.3) 

The initial state of the apparatus and the bath is assumed to be factorized at some remote initial time t^. At that time 
the bath was in equilibrium at temperature /?~^ while the apparatus was in an arbitrary state. Explicit integration 
of Eq. (D.3), through 

(t) = 6n (io) e-'"'"(*-*«) - ic,;, /* dt' a (f) e'f^'"(*-*'), (D.4) 



J to 

allows us to eliminate the bath and to write a closed, exact equation of motion for the operator a. In the weak-coupling 
regime it is convenient to go the rotating frame by means of the transformation 

a (0 = a (0 e'"(*-*«^ (D.5) 

From Eqs. (D.2), (D.4) and (D.5) we obtain 

^ = - E r « (0 e'("-"'")(*-*') +C{t), (D.6) 

where C (t) appears as a gaussian quantum noise defined by 

C it) = -lY^CmU (to) e'("-^-)(*-*«). (D.7) 
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This noise is characterized by the properties 



[CW,C^(i')] = El^l'^''""''"'^^*"*'^' [C(i).C(i')] = 0, (D.8) 



2 .i(Q_0„)(t-t') 



(CHf) C m = E '"";,na^_i ' (C (*) C it')) = 0. (D.9) 

m 

We are interested in the thermodynamic limit for the bath, where the number of bath modes goes to infinity while 
their frequencies 0^ tend to a continuum. For simplicity, we assume that the values Clm have a constant spacing 
A — > 0, and that the coupling \cm\ is a vanishingly small constant: 

(D.IO) 

V 71 

We then have 

E e'("-"-H*-*') = 2jS {t - t') . (D.ll) 

m 

(If the spectrum of Q,„i is bounded, the (5-fimction stands for a narrow peak, with a width smaller than all the 
characteristic times of the problem if the range of is sufficient.) Eq. (D.ll) ensures that the retardation effect 
can be neglected in Eq. (D.6), which becomes: 

^ = -ra(t) + at), (D.12) 
and that the properties of the noise are simplified into 

[C(t),CUO] = 27-5 (*-*'), [CW,C(*')]=0 (D.13) 



(C^ (0 c w) = lj X-^-V <^ ^ = °- 

Apart from the very short time-scale involved in the noise ^ {t) , two time-scales, and 7""'^, enter the dynamical 
equations (D.5), (D.12) for the operator a{t) generated by the apparatus-bath Hamiltonian Ha,b- Other, longer 
time-scales will also be induced by the interaction with the tested system. We assume that the dynamical frequency a 
is the largest characteristic frequency of the problem. The transformation (D.5) then acounts for the fast, oscillatory 
motion of a{i). The evolution of a (t) takes place on larger time-scales, of order 7"^ or more. We can thus expand 
Eq. (D.14) for large a according to 



(Ct(OC(i)) «^/da;e'"(*-Og 

•' n=0 



da J e/5^« - 1 

(D.15) 



„ n! V da 

Since this expression will be integrated over a function of time which varies on the scale 7, its successive terms yield 
an expansion in powers of j/a, and we can thus replace Eq.(D.14) in the limit a 7 by 

<C+ it') C it)) = 27^1^, (C it) C it')) = 0. (D.17) 

If the noise is redefined as 

6o(i) = 4=CWe-^"(*-*°\ (D.18) 
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Eq. (D.12) reduces to the equation of motion (2.21) for oq (t) , and Eqs. (D.13), (D.16) coincide with the conditions 
(2.23), (2.24), (2.25) on bg (t) . We can hkewise recover Eqs. (2.21-2.25) for i and fc ^ by adding Wfc to a and by 
couphng all the modes i of the apparatus with the modes m of the bath. 

If we had not used the approximation a ^ 7, we would have obtained the evolution of the average particle number 



no (i) = (a^ (t) a (i)) in the lowest level of the apparatus by integration of Eqs. (D.5), (D.12), using (D.14) instead of 



(D.17). This yields 



no (t) = no (to) e-^-(*-o) + 1 f % 1 - ^^^^-^^)^^-^^) " (d.19) 



TT 



(w — a) + 7^ 



2 1 



instead of Eq. (2.27). In particular, we find for arbitrary 7/a that uq (t) relaxes for large times towards 



no {t)^l I (D.20) 



It is only for 7 ^ a that Eq. (D.l£) reduces to Eq. (2.27), and that no {t) reaches the expected equilibrium Bose 
factor 1/ (e^''" — l) for large times. 

Altogether the weak-coupling condition 7 <C a, which will be enforced throughout this paper, ensures that the 
quantum Langevin equation (2.21) can be obtained from the apparatus-bath Hamiltonian (D.l). Owing to this 
condition, while the coupling parameter 7 governs the dynamics of the apparatus during its relaxation to equilibrium, 
the equilibrium properties of the apparatus remain unaffected by the presence of the bath. Note finally that the 
simplifying assumptions we made on the bath (equally spaced levels constant |c„i|) become irrelevant in the 
weak-coupling limit. 
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